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An investigation is performed of the Lorentz- violating electrodynamics extracted from the renor- 
malizable sector of the general Lorentz- and CPT-violating standard-model extension. Among the 
unconventional properties of radiation arising from Lorentz violation is birefringence of the vacuum. 
Limits on the dispersion of light produced by galactic and extragalactic objects provide bounds of 
3 x 10~ 16 on certain coefficients for Lorentz violation in the photon sector. The comparative spectral 
polarimetry of light from cosmologically distant sources yields stringent constraints of 2 x 10~ 32 . 
All remaining coefficients in the photon sector are measurable in high-sensitivity tests involving 
cavity-stabilized oscillators. Experimental configurations in Earth- and space-based laboratories are 
considered that involve optical or microwave cavities and that could be implemented using existing 
technology. 



I. INTRODUCTION 

Lorentz symmetry underlies the theory of relativity 
and all accepted theoretical descriptions of nature at the 
fundamental level. A crucial role in establishing both 
the rotation and boost components of Lorentz symmetry 
has been played by experimental studies of the proper- 
ties of light. In the classic tests, rotation invariance is 
investigated in Michelson-Morley experiments searching 
for anisotropy in the speed of light, while boost invariance 
is studied via Kennedy-Thorndikc experiments seeking a 
variation of the speed of light with the laboratory velocity 

In this work, a theoretical study is performed of vari- 
ous experiments testing Lorentz symmetry with light and 
other electromagnetic radiation. The analysis is within 
the context of the Lorentz- and CPT-violating standard- 
model extension Q, developed to allow for small gen- 
eral violations in Lorentz and CPT invariance J5|. The 
lagrangian of this theory includes all observer Lorentz 
scalars formed by combining standard-model fields with 
coupling coefficients having Lorentz indices. At the level 
of quantum field theory, the violations can be regarded as 
remnants of Planck-scale physics appearing at attainable 
energy scales. The coefficients for Lorentz violation may 
be related to expectation values of Lorentz tensors or 
vectors in an underlying theory ||. To date, experimen- 
tal tests of the standard-model extension have been per- 
formed with hadrons if ^Lol, protons and neutrons [fj"T|| , 
electrons photons |0 @ , and muons Jl|]. 

In the present context of studies of electrodynamics, 
the standard-model extension is of interest because it 
provides a general field-theoretic framework for inves- 
tigating the Lorentz properties of light. The theory 
contains as a subset a general Lorentz-violating quan- 
tum electrodynamics (QED), which includes a general 
Lorentz-violating extension of the Maxwell equations. 
We study experiments that can measure the coefficients 



for Lorentz violation in this generalized electrodynamics. 
Our attention is restricted here to exceptionally sensitive 
experiments that could be in a position to detect the mi- 
nuscule effects motivating the standard-model extension. 

A basic feature of Lorentz-violating electrodynamics is 
the birefringence of light propagating in vacuo. This re- 
sults in several potentially observable effects, including 
pulse dispersion and polarization changes. One goal of 
this work is to consider the implications of these effects 
for the propagation of radiation on astrophysical scales. 
We use available observations to constrain certain coeffi- 
cients for Lorentz violation. 

Another goal of this work is to analyze modern ver- 
sions of some classic tests of special relativity based on 
resonant-cavity oscillators which have extreme 

sensitivity to the properties of electromagnetic fields. 
These experiments depend on the Earth's sidereal and 
orbital motion. However, the advent of the International 
Space Station (ISS) makes it feasible to perform labora- 
tory experiments in space, where the orbital motion can 
yield different sensitivity to Lorentz-violating effects [^0) . 
We consider here both space- and Earth-based laboratory 
experiments with resonant cavities. 

The structural outline of this paper is as follows. Sec- 
tion [n] presents some basic results and definitions for 
the general Lorentz-violating electrodynamics and out- 
lines the connection to some test models. We then con- 



sider birefringence experiments, beginning in Sec. [II A 
with some general issues. Constraints stemming from 
the resulting effects on pulse dispersion from astrophysi- 
cal sources are addressed in Sec. IIIB, while those from 



polarization changes over cosmological scales are treated 
in Sec. 



IIIC 



A general analysis for laboratory-based 
exp erimen ts on the Eart h a nd in space is presented in 
Sec. IV A . Sections [VB| and |IV C apply this analysis to 
experiments with optical and microwave resonant cavi- 
ties. We summarize in Sec. [v|. Throughout this work, 
we adopt the conventions of Ref . m . 



II. LORENTZ- VIOLATING ELECTRODYNAMICS 

This section provides some background and contextual 
information about the general Lorentz-violating electro- 
dynamics. The basic formalism is presented, and some 
definitions used in later sections are introduced. We also 
discuss the connection between this theory and some test 
models for Lorentz violation. 



A. Basic Theory 

The standard model of particle physics is believed to 
be the low-energy limit of a fundamental theory that 
includes all the forces in nature. The natural scale of 
this fundamental theory is likely to be determined by the 
Planck mass. The possibility that Lorentz- and CPT- 
violating signals from this theory may be observable at 
energies attainable today led to the development of the 
standard-model extension Q, which is a general theory 
based on the standard model but allowing for violations 
of Lorentz and CPT symmetry Q . The additional terms 
must be small because the usual standard model agrees 
well with experiment. They may originate from spon- 
taneous symmetry breaking in the fundamental theory 

The standard-model extension can be defined as the 
usual standard-model lagrangian plus all possible ad- 
ditional Lorentz- and CPT-violating terms involving 
standard-model fields that maintain invariance under 
Lorentz transformations of the observer's inertial frame. 
This invariance ensures that the physics is independent of 
the choice of coordinates. The Lorentz violation is asso- 
ciated with rotations and boosts of particles or localized 
field configurations in a fixed observer inertial frame. 

Many of the detailed investigations of the standard- 
model extension have been performed under the simpli- 
fying assumption that the additional Lorentz- and CPT- 
violating terms preserve the SU(3)xSU(2)xU(l) local 
gauge symmetry of the usual standard model. Another 
widely adopted simplifying assumption is that the coeffi- 
cients for Lorentz violation are independent of position. 
This implies the violation is restricted to the Lorentz 
symmetry instead of the full Poincare symmetry and has 
several useful consequences for experiment, including the 
conservation of energy and momentum. It is also of- 
ten convenient to restrict attention to the renormalizablc 
sector of the theory, since this is expected to dominate 
the physics at low energies. However, nonrenormalizable 
terms are known to play an important role at higher en- 



renormalizable terms, which involve operators of mass 
dimension four or less. The relevant lagrangian is H 



ergies 21 



Extracting terms involving the photon fields from the 
standard-model extension yields a Lorentz- and CPT- 
violating extension of QED Q . The fermion sector of this 
theory has been widely studied. Here, we focus atten- 
tion on the pure-photon sector and limit attention to the 



C = -\F^F^ + \{k AF Ye KX ^A x F^ 



4 \ K F ) nXfiv 



(1) 



where = d^A v — d v A^. This theory maintains the 
usual U(l) gauge invariance under the transformations 
qA^ — > qA^ + d^A. The lagrangian contains the stan- 
dard Maxwell term and two additional Lorentz-violating 
terms. The first of these extra terms is CPT odd, and its 
coefficient (kAF) K has dimensions of mass. The other is 
CPT even. Its coefficient {kF)n,\p,u is dimensionless and 
has the symmetries of the Riemann tensor and a vanish- 
ing double trace, which implies a total of 19 independent 
components. 

The CPT-odd term has received much attention in the 
literature p2| . This term provides negative contributions 
to the canonical energy and therefore is a potential source 
of instability. One solution is to set the coefficient to zero, 
(fc^f) K = 0. This is theoretically consistent with radia- 
tive corrections in the standard-model extension and is 
well supported experimentally: stringent constraints on 
Uaf have been set by studying the polarization of radia- 
tion from distant radio galaxies . 

In contrast, much less is known about the CPT-even 
coefficient hp. Theoretical studies show that it provides 
positive contributions to the canonical energy and that 
it is radiatively induced from the fermion sector in the 
standard-model extension P,p3[. Constraints on some 
components have recently been obtained from optical 
spectropolarimetry of cosmologically distant sources |l5| . 
In the present work, we focus on the experimental impli- 
cations of this CPT-even term. The coefficient (A;af) k is 
set to zero for the analysis. 

The equations of motion from lagrangian (|l]) are 



d a F° + (M MQ/37 d Q i^ = o. 



(2) 



These are modified source-free inhomogeneous Maxwell 
equations. The homogeneous Maxwell equations, 



d,F^ = | 



e^d^x = 0, 



(3) 



remain unchanged. 

Although it lies beyond our present scope, the tech- 
niques presented here and the results obtained can be 
generalized to the nonrenormalizable sector. The non- 
renormalizable terms can be classified according to their 
mass dimension. The dimensions of the corresponding 
coefficients are inverse powers of mass, and it is plausi- 
ble that these coefficients are suppressed by correspond- 
ing powers of the Planck scale. Terms of this type ap- 
pear in various special Lorentz-violating theories, includ- 
ing noncommutative field theories incorporating QED 
PjJ. Indeed, any coordinate- independent theory with 



a photon sector containing nonrenormalizable Lorentz- 
violating terms must be a subset of the standard-model 
extension. It would be interesting to provide a detailed 
study of the nonrenormalizable terms in the Lorentz- 
violating electrodynamics and their experimental signals. 



B. Analogy and Definitions 

A useful analogy exists between the Lorentz-violating 
electrodynamics in vacuo and the conventional situation 
in homogeneous anisotropic media 0] . The idea is to de- 
fine fields D and H by the six-dimensional matrix equa- 
tion 



1 + k D e K D B 
k H e 1 + khb 



(4) 



where E and B are the electric and magnetic fields ob- 
tained from solving the modified Maxwell equations (|^). 
The 3x3 matrices kde , K h b, kdb, and kh e are defined 

by 





-2(fc F ) 0j0fe , 




{K HB y k = 


±e jpq e krs (k F ) pqrs , 




(K DB ) jk = 


-(K HE ) kj = (k F f ]pq e kpq . 


(5) 



The double-trace condition on (k F ) translates to 
the tracelessness of (kde + khb), while (k F ) K [\ tlv ] = 
implies the tracelessness oi kdb = — (khe) T - This 
leaves kde and khb with eleven independent elements 
and the matrix kdb — — (khe) T with eight, which to- 
gether represent the 19 independent components of k F . 
Note also that kde and khb are parity even, while 
kdb — —{ k he) T is parity odd. 

With these definitions, the modified Maxwell equations 
(||), (||) take the familiar form 



V x H - d D = 0, V • D = 0, 

V x E + d B = 0, V • B = 0. 



(6) 



As a consequence, many results from conventional elec- 
trodynamics in anisotropic media also hold for this 
Lorentz-violating theory. For example, the energy- 
momentum tensor takes the standard form in terms of E, 
B, D and H. This implies the usual Poynting theorem, 
which can be applied in conjunction with the symme- 
tries of the matrices in Eq. (||) to show that the vacuum 
is lossless. 

For the applications to be addressed in later sections, 
it is convenient to introduce the following decomposition 
of (fcF)fcA/ii/ coefficients: 



{k e+ y k = \(k D e + K H B) 

(k e _y k = \(k de - K HB y k - \V k {KDE) U , 



(k 0+ y k = \{k db + K HE y k , 
(k -) jk = \{k DB - K HE y k i 

ktr = \ {kde) U ■ 



(7) 



The first four of these equations define traceless 3x3 
matrices, while the last defines a single coefficient. All 
parity-even coefficients are contained in k e+ , K e _ and K tT , 
while all parity-odd coefficients are in k Q+ and k q _ . The 
matrix k 0+ is antisymmetric while the other three are 
symmetric. 

The form of this decomposition helps in determin- 
ing the portion of the parameter space to which ex- 
periments are sensitive and how different experiments 
might overlap. For example, typical laboratory experi- 
ments with electromagnetic cavities search for rotation- 
violating parity-even observables. The sensitivity of such 
experiments is therefore expected to be dominantly to the 
ten rotation-violating parity-even coefficients k e+ and 
k e -. For those observables depending at leading order 
on the velocity, the eight coefficients k a+ and k q _ can be 
expected to play a role. Finally, at second order in the 
velocity one can expect the sole rotation-invariant quan- 
tity K tr to affect measurements. These considerations are 
confirmed by the results of the detailed analysis in the 
sections below. 

As another example of the use of the decomposition 
(0), recall that birefringence is known to depend on ten 
linearly independent combinations of the components of 
k F , which can be chosen as |13] 



((fcF)° 213 , 


(fcF)° 123 , 






(fc F ) 0202 - 


(fcF) 1313 , 


(k F ) 0303 


" (fcF) 1212 


(fcF) 0102 4 


(fcF) 1323 , 


(M 0103 


- (fc F ) 1223 


(k F ) 0203 + 


(fcF) 12 ", 


(fcF) " 2 


+ (fcF)° 323 


(fcF) " 3 " 


(fcF)° 223 , 


(fc F )° 212 


- (fcF)° 313 



Relating these to the k matrices, we find 



(£ e +) 



jk 



(Ho-) 



jk 



-(k 3 + k 4 ) k 5 k 6 
k 5 k 3 k 7 

k 6 k 7 fc 4 

2fc 2 —k 9 fc 8 
~k 9 -2k 1 k w 
k s k w 2(fc 1 - fc 2 ) 



(9) 



In this way, we can see directly that birefringence is con- 
trolled by the matrices k e+ and k _. 

In terms of the k matrices defined in Eq. (|[), and as- 
suming as before that (^af)" = 0, the lagrangian (|l|) 
becomes 

C = \(E 2 - B 2 ) + \E ■ (k DE ) -E-\B- (khb) • B 



+E ■ (kdb) ■ B. 



(10) 



Similarly, using instead the k matrices defined in Eq. (Q), 
we find 



£ = ±[(l + k tI )E 2 -(l-k tr )B 2 } 



\E ■ {k e 



E 



B. 



B 

(11) 



The form of Eq. (|TT]) shows that a nonzero coefficient 
Ktr shifts the effective permittivity e and effective perme- 
ability /j, by (e — 1) = — (/J -1 — 1) = Ktrj corresponding 
to a shift in the speed of light. However, it is possible 
to remove an overall shift in the speed of light by mak- 
ing advantageous coordinate transformations accompa- 
nied by suitable field redefinitions, which combine to set 
e = = 1 and transfer the Lorentz violation to a dif- 
ferent sector of the theory. An explicit example of this 
procedure is provided in the next subsection for a toy 
model involving scalar QED. 

In the general context of the standard- model extension, 
such transformations modify various other coefficients for 
Lorentz violation. In fact, similar transformations can 
move the nine independent coefficients k e ^, and k tI 
into other sectors of the theory. Note that this effect 
is frame dependent because the coefficients mix under 
boosts. Note also that the possibility of absorbing K e _, 
k +, Ktr elsewhere offers insight as to why birefringence 
experiments, which directly compare light with light, are 
insensitive to these coefficients. However, cavity experi- 
ments involve comparisons of radiation with matter, so 
all 19 coefficients are observables in this case. 



C. Connection to Some Test Models 

Several phcnomcnological test models for Lorentz 
properties of light have been proposed. The standard- 
model extension contains all observer-independent 
sources of Lorentz violation in terms of known particles, 
so it is expected to incorporate the existing test models 
as special cases. In this subsection, we comment on the 
relationships to some popular test models. 

Since typical test models assume only one type of mat- 
ter other than the photon, it suffices for our purposes to 
consider a toy version of the standard-model extension 
that includes only one scalar field and a limited type of 
Lorentz violation. We therefore work with a model of 
Lorentz-violating scalar QED, defined by the lagrangian 



C = {if + (k+fXD^Dut - m 2 
" \{k F ) KX ^F KX F^. 



(12) 



In this expression, the covariant derivative takes the 
usual form, D^<f) — d^c/) + iqA^cj), and for simplicity we 
have limited the types of Lorentz violation to those de- 
scribed by a real symmetric coefficient {k^)^ u and by a 
coefficient (kp)^^ of the type in Eq. ([!]). 

An interesting test model for Lorentz violation is pro- 
vided by the kinematical framework of Robertson ]2q| 



and its extension to arbitrary synchronizations by Man- 
souri and Sexl Q. These approaches suppose the ex- 
istence of a "preferred" frame in which light propagates 
isotropically as measured by a standard set of rods and 
clocks. The Lorentz transformation between observers is 
then generalized to incorporate small changes from the 
conventional boosts in special relativity. Within a given 
synchronization, three parameters g$, g\, gi are needed 
to fix the generalized Lorentz transformation and hence 
to characterize the Lorentz violation. 

The construction of the generalized Lorentz transfor- 
mation can be illustrated in the context of the model 
(|l2l). Consider the special case of the model for which 
only the coefficient (fc^) 00 is nonzero in a certain frame 
S. Writing this coefficient as (fc^) 00 = k 2 — 1, where k 2 
deviates slightly from 1, the lagrangian takes the form 



C 



+ UE 2 



^D^<j)+{k 2 
B 2 ). 



l)|A)C 



m 2 (jy <p 



(13) 



In the £ frame, the propagation of light is rectilinear and 
isotropic, so it may be identified with the preferred frame 
of the test model. The Lorentz violation appears only in 
the <p sector of the lagrangian, which we can suppose de- 
scribes the detailed physics of the rods or clocks in the 
test model. 

The generalized Lorentz transformations T^ v consid- 
ered in the kinematical test models are the linear trans- 
formations x'^ = v x v from the preferred frame S to 
a coordinate system S attached to an observer moving 
at constant velocity in the preferred frame. By construc- 
tion, the observer S defines coordinates using the same 
rods and clocks and a prescribed synchronization. How- 
ever, in the present context the Lorentz-violating prop- 
erties of the rods and clocks are fixed by the Lorentz- 
violating scalar term in the lagrangian dH). The general- 
ized Lorentz transformations T^ v from £ to S are there- 
fore also determined in the context of the model jl^). 
They are the transformations leaving invariant the scalar 
sector and hence preserving the combination rf v + (fc^)^ 
up to a possible resynchronization. For example, for 
the special case of Eq. (p"3|), the Robertson parameters 
are found to be g = l/g 1 = - /3 2 )/(l - fc 2 /3 2 ), 
92 = I- The corresponding Mansouri-Sexl parame- 
ters arc a = 1/6 = — k 2 /3 2 ), d = 1, with e = 
— /3(1 — k 2 (3 2 )/(l — (3 2 ) in Einstein synchronization or 
e = —k 2 (3 in slow-clock synchronization. In contrast, the 
standard Lorentz transformations preserve rf v ' . 

In this simple example, the transformation leaves 
invariant the rods and clocks, while A^„ leaves invari- 
ant the speed of light. Both are equally valid. In the 
frames related by T^ v , observers agree on rod lengths and 
clock rates but disagree on the velocity of light. More- 
over, the velocity of light is no longer isotropic as mea- 
sured by these rods and clocks. In contrast, observers 
related by Lorentz transformations agree that light prop- 



agates isotropically with speed 1 but may disagree on rod 
lengths and clock rates. The description is a matter of 
coordinate choice, and one can move freely from one to 
the other using T M „, A^ v , and their inverses. 

Note that a "preferred" frame in which light prop- 
agates isotropically typically fails to exist in the full 
standard-model extension, although in principle one can 
impose the existence of such a frame by suitably re- 
stricting the coefficients for Lorentz violation. From this 
perspective, the special status enjoyed by photons rel- 
ative to other particles in the kinematical test models 
appears somewhat unnatural, and the structure of the 
standard-model extension offers more general possibil- 
ities for kinematical frameworks. Note also that the 
standard-model extension addresses modifications to all 
known particles, so the effects on physical rods and clocks 
can be directly analyzed. This is infeasible in kinematical 
frameworks, which consider the transformations between 
frames rather than the underlying physics. 



Another interesting test model is the c 2 model 27 1, de- 
veloped for application to studies of Lorentz invariancc 
as a limiting case of the THep formalism [ p8|]29| . The c 2 
model is defined by a lagrangian describing the behav- 
ior of classical pointlike test particles in the presence of 
electromagnetic fields. The model assumes the existence 
of a "preferred" frame in which the limiting speed of the 
test particles is 1, while the speed of light is c. 

To see the relation between the c 2 model and the model 
( |l2| ) , consider another lagrangian written in a frame S as 

C = {D^D^cj) - mV0 + \[E 2 - k 2 B 2 ), (14) 

where k 2 deviates slightly from 1 as before. In this theory, 
the Lorentz violation appears in the photon sector. With 
the identification k — c, the lagrangian for this sector is 
identical to that of the c 2 model. Moreover, the </> sector 
is conventional, representing a quantum field theory of 
minimally coupled scalar particles. The model (|l4|) can 
therefore be regarded as the field-theoretic equivalent of 
the c 2 model. 

The two models ([TI]) and ( |l3| ) are related by the coordi- 
nate transformation t — > t/k, x — > x followed by the field 
redefinition — > A^/k and charge rescaling q — > kq. 
They therefore describe the same physics. Although it 
is possible to choose coordinates so that either the pho- 
ton or the scalar propagates conventionally, the Lorentz 
violation cannot be eliminated simultaneously from both 
sectors. 

We thus see that the c 2 model is contained in the the- 
ory ( p^ ) as a special case. In the terminology of Eq. ([?]), 
the parameter c 2 could be identified with the combina- 
tion of coefficients (1 — k tr )/(l + kti), as can be seen 
from Eq. (|ll|). However, caution is required in interpret- 
ing bounds obtained with the c 2 model in terms of Rt r 
because the identification is valid only in a frame S with 
conventional particles, which typically fails to exist in the 
standard-model extension. 



III. ASTROPHYSICAL TESTS 

In this section, we consider radiation propagating in 
free space. The Lorentz-violating electrodynamics pre- 
dicts birefringence, which allows sensitive tests of Lorentz 
symmetry from observations of radiation propagated over 
astrophysical distances. We begin with some general the- 
ory, and then we obtain two sets of bounds on Lorentz 
violation from velocity and birefringence constraints. 



A. General Theory 

The basic features of plane-wave solutions to the 
Lorentz-violating electrodynamics have been presented 
in Refs. |4|,[l5|], so only relevant essentials are given here. 
With the standard ansatz F fiv (x) — F^ u (jj)e~ WaX for a 
plane wave with wave 4- vector p a = (p°,p), the equation 
determining the dispersion relation and the electric field 
E is the modified Ampere law 



M jk E k _ ^s:> k p 2 -p> P k - 2{k F y^ k pp Pl )E k = 



(15) 



The dispersion relation is obtained as usual by requiring 
vanishing determinant of M^ k . It suffices for our present 
purposes to consider only leading-order effects in the co- 
efficients (kp) K \^ for Lorentz violation. To leading or- 
der, one finds 



(l+p±a) 



where 



with 



w*f> = (k F ) a ^p^, P »=p»/\p\. 



(16) 



(17) 



(18) 



Note that p^p and fPa are observer Lorentz scalars, which 
implies p and a are scalars under observer rotations. 

The dispersion relation ( |l6| ) has two solutions, with 
corresponding electric fields E± . In conventional electro- 
dynamics, the dispersion relation is p° — \p\ and all fields 
E perpendicular to p are solutions, so the propagation is 
independent of the polarization. However, in the present 
case the propagation is goverened by two specific modes 
E±, with the general solution to Eq. ( p^ ) being any lin- 
ear combination of the two. This leads to birefringence: 
light generically has two components, each propagating 
independently. 

There are several possible definitions for the velocity of 
the radiation, including the phase velocity v J p = p^pi jp 1 , 



the group velocity v g = (V^p , and the velocity of en- 



on 



ergy transport v\ = 9 J{ '/9 W , where is the energy- 
momentum tensor. With the analogy discussed in Sec. 



II B, one can show by standard arguments that v g = v e 



for a wave with fixed p. Also, Eq. (16) can be used 
to find explicit leading-order expressions for the mag- 
nitudes of the phase and group velocities. We thereby 
obtain v p = v g = v e = 1 + p ± a to leading order in 
(fc_F)«A^i/- Note also that, to leading order, we can write 
p a = (p°,p)/\p\ ~ hi the expressions ( |l7|) for p 

and a. The quantity v can be regarded as the direc- 
tion of propagation of the radiation, since the difference 
between it and the other velocities arises only at higher 
order and is irrelevant here. 

The mode dependence of the velocity offers interest- 
ing possibilities for experimental tests of the theory. The 
velocity difference is 



Av = v+ — V- = 2cr, 



(19) 



and is expected to be tiny. However, for sufficiently large 
path lengths this difference might become apparent in 
the form of observable effects on the pulse shape or the 
polarization of radiation. In the next two subsections, we 
exploit these features to obtain constraints on {kp) K \^ v . 

An explicit form for the solutions E± is needed for 
some of the analysis. Using the dispersion relation (16), 
the matrix in the Ampere law (Q5h can be written 



Mf = 



■[2(p±a)S jk 



- p>p k 



2F fe ' 



V 



(20) 



The form of this matrix shows that the solutions E± are 
wavelength independent but vary with the direction of 
propagation. Also, at leading order in {kF) K \^v the dif- 
ference M + — M_ is proportional to the identity, so the 
leading-order solutions E+ and E- are perpendicular. In 
fact, at leading order, E± are perpendicular to p as well. 

To express E± explicitly, a choice of inertial frame must 
be made. It is convenient to adopt a standard reference 
frame to report the results of observations and hence ul- 
timately to place constraints on the set of coefficients Uf- 

A natural choice for the reference frame is a Sun- 
centered celestial equatorial frame with the Z-axis 
aligned along the celestial north pole at equinox 2000.0. 
The Z-axis is then at a declination of 90°, and the X- 
and y~axes lie at declination 0° and can be chosen to 
be at right ascension 0° and 90°, respectively. The unit 
vector X thus points towards the vernal equinox on the 
celestial sphere. The time T is chosen such that T = 
when the Earth crosses the XY plane on a descending 
trajectory. In what follows, we adopt this standard frame 
to report results. 

For the practical determination of E± for a given wave, 
it is easiest first to work in a special 'primed' frame chosen 
for that wave. The result of the calculation can then be 
related to the standard Sun-centered frame by perform- 
ing a suitable observer rotation. A convenient primed 



frame for a given wave is the frame in which the wave 
4- vector takes the form p la — (1; 0, 0, 1) to leading order. 
The solution for E± can be expressed explicitly in terms 
of the coefficients k' F in this frame. Up to a normaliza- 
tion, it is found to be E± oc (sin£,±l — cos£,0), where 
tan£ = 2k' 12 /{k' 11 - k' 22 ). The two modes are thus 
linearly polarized. 

From the solutions E± and dispersion relation, it is 
evident that a and £ are the relevant parameters for 
birefringent effects for a particular source. In particu- 
lar, crsin^ ~ k' 12 and a cos £ = ^{k' 11 — k' 22 ) represent 
the minimal linear combinations of k' F that govern bire- 
fringence. The parameter p is common to both modes, 
but does not contribute to birefringence and cannot be 
detected in the experiments discussed below. 

The results in the primed frame can be related to the 
standard frame by a suitable observer rotation, described 
in Appendix [A]. The direction of travel of the light in 
the standard frame determines two vectors ^ and c? in 



k a space (see Eq. (A4)), and it turns out that the bire- 
fringence of the light depends on the two specific linear 
combinations of the coefficients k a in Eq. (||) that are 
parallel to these vectors. 



B. Velocity constraints 

For the two radiation modes E± propagating over a 
distance L, the velocity difference (|l9|) induces a differ- 
ence At ss AvL between the two travel times. Local mea- 
surements made on radiation emitted as a single burst 
from a distant source can therefore provide sensitivity to 
the coefficients k a for Lorentz violation p0| . 

To apply this idea, it is useful to consider distant 
sources that produce radiation in a relatively narrow 
burst characterized by a small width w, such as mil- 
lisecond pulsars or sources of gamma-ray bursts. These 
sources typically produce essentially unpolarized radia- 
tion, so the intensity of each mode should be compara- 
ble. The burst can then be regarded as a superposition 
of two independently propagating pulses, one for each 
mode. For a sufficiently great distance L, a nonzero Av 
would cause the two pulses to separate enough to become 
distinguishable. This type of signal would manifest itself 
as two pulses with similar time structure but differing in 
arrival time. The pulses would each be linearly polarized, 
and they would have mutually perpendicular polarization 
angles. 

If only a single pulse is observed, a limit on Lorentz 
violation can be deduced. The relationship between the 
observed pulse width w and the source pulse width w s 
is approximately w Q ~ w s + At. Observations of w 
can therefore be used to obtain a conservative bound on 
At = AvL = 2aL and hence a bound on the coefficients 
k a . 
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Table 1. Source data for velocity constraints. 

Table 1 lists data for fifteen sources suitable for placing 
this type of constraint. The first five lines list gamma- 
ray bursts with known redshifts. The widths listed for 
these contain all significant time structure of the pulse. 
The distance L is determined from the redshift by the 
look-back time in a conservative cosmology for a matter- 
dominated universe with Hubble constant Hq = 80 km 
s _1 Mpc -1 . The next eight sources in the table are 
millisecond pulsars. The listed pulse width is that at 
10% peak intensity. The final two sources are giant-pulse 
pulsars. These exhibit intense pulses with characteristic 
widths on the order of several ^s. 

For each source in Table 1, we take a < w /2L as a 
bound on a. For a single source, this places constraints 
on a two-dimensional subspacc of the full 10-dimcnsional 
parameter space of the coefficients k a . The subspace is 
represented by the linear combinations a sin £ and a cos £ 
associated with that particular source. To bound all ten 
coefficients fc a , ten linearly independent constraints of 
this type are needed. This is feasible using five or more 
sources at different positions on the sky. 

We proceed by assuming the constraint for each source 
in Table 1 is consistent with a measurement of a = 0, and 
we take the bound a < w /2L as a reasonable estimate 
of the error in a null measurement. The associated % 2 
distribution is \ 2 = ^L 2 a 2 /w 2 , where the sum is over 
the fifteen sources. This is a quadratic form in k a . Con- 
sidering |fc a | and minimizing \ 2 with respect to the other 
nine degrees of freedom, we obtain a bound of 

-16 



\k a \ < 3 x 10" 



(21) 



in the Sun-centered celestial equatorial frame, at the 90% 
confidence level. 

This bound is much less stringent than that obtained 
through polarization measurements, as discussed below. 
However, the method is relatively straightforward and 
avoids some of the complexities involved in the polariza- 
tion analysis. 



C. Polarization constraints 



In this subsection, we expand on the material found 
in Ref. p"5| . An improvement on the previous result is 
made by considering the cosmological redshift of light. 

A general electric field E can be decomposed into 
its birefringent components E±. Defining unit vectors 
e± = E±/\E±\, the decomposition is 



E(x) = {E+e+e- 1 ^ + E-e-e-^-ty 



(22) 



The differing phase velocities of the two modes results in 
a change in relative phase as the wave propagates, given 

by@ 

A(f> = (p° + - p°_)t ps 2irAv p L/\ ps 47rj£/A, (23) 

where Av p is the difference in phase velocities, A is the 
wavelength, and L is the distance traveled. The phase 
change modifies the polarization state of the radiation, 
with larger effect for more distant sources. Appendix [b] 
provides a brief review of pertinent concepts involving 
polarization in the present context. 
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FIG. 1. Rotation of the Stokes vector about sj 



In the primed frame described in Sec. Ill A, the Stokes 
vectors for e± are s± — ±(cos£,sin£,0). These vec- 
tors correspond to opposite points on the equator of the 
Poincarc sphere, as expected for linearly polarized modes. 
As described in Appendix [b], the axis of rotation induced 
by the phase change A0 = 4iraL/\ is therefore in the s 1 - 
s 2 plane. This affects both %j) and x> as can be seen from 
Fig. @. 

A change in phase can arise from a change in either 
L or A. The induced change in the polarization depends 
not only on hp but also on the initial polarization. For 



cosmological sources, it may be impossible to determine 
independently the polarization at the source, in which 
case one cannot determine whether a change in polariza- 
tion is strictly due to a change in L. It is therefore of 
more interest to focus on the wavelength dependence of 
the polarization change. Making the reasonable assump- 
tion that the emitted polarization is relatively constant 
over a given range of wavelengths, the relevant quantity 
becomes the phase shift as a function of wavelength, 



(24) 



relative to a reference wavelength Ao- Standard spec- 
tropolarimetric techniques then allow a measurement of 
this effect. Note that knowledge of physical processes in 
certain classes of objects producing the polarized radi- 
ation might make it feasible to include a known initial 
polarization in the analysis, but this is unlikely to im- 
prove significantly the constraint obtained here. 

The effect on the measured polarization as the wave- 
length is changed can be visualized using the Poincare 
sphere. Suppose a source produces radiation with con- 
stant polarization over a range of wavelengths. This radi- 
ation can be represented by a single point on the Poincare 
sphere. As the light propagates towards the Earth, the 
presence of Lorentz violation causes this point to rotate 
along an arc on the sphere. For any fixed wavelength, the 
rotation axis and rate depend on the coefficients (kp)^^ 
and on the position of the source on the sky. However, 
Eq. (M) shows that shorter wavelengths rotate more than 
longer ones. Therefore, as measurements of the Stokes 
vector are made over a range of wavelengths, the results 
trace a circular arc on the surface of the Poincare sphere. 

Let ipo and \o represent the observed polarization of a 
point on this arc with the reference wavelength Ao. Using 
this point as a reference, we require the change in polar- 
ization relative to this point induced by Eq. (fi4|). This 
polarization change is given by 



=m^(<^)s fc (</>o,Xo), 



(25) 



where m? is the rotation matrix about s + by S(f>. The 
matrix m 3 is analogous to the Mueller matrix used in 
polarimetry to describe the effects of various filters and 
polarizers on light. Its explicit form is given as Eq. (8) 
of Ref. p5| . The angle x> which controls the amount 
of circular polarization, is absent from most published 
spectropolarimetric data. It is therefore most effective to 
focus attention on the change 5ip — ip — tJjq m ip from the 
reference value i/>o, which is given as Eq. (9) of Ref. [JL5). 
A procedure for fitting this equation to existing spec- 
tropolarimetric data is also provided in this reference, 
and a 90% confidence-level bound of \k a \ < 3 x 10~ 32 is 
obtained from spectropolarimetric data for 16 cosmolog- 
ical sources. 

In the present work, we use the same procedure to 
obtain a slight improvement on the existing bound by 



incorporating the redshift of the light as it propagates 
to the Earth. Cosmological redshift implies that over 
the path traveled the light has shorter wavelength than 
observed. Taking the same conservative cosmology as in 
the previous subsection and integrating the phase change 
over the propagation time yields 



Acj) = 



47TCT 2 



1 



1 



(26) 



where A is the observed wavelength and z is the redshift. 
To account for the redshift, it therefore suffices to replace 
L with L c ff = 2(1 — (1 + z)~ x / 2 )/ H in the analysis. 



Source 



IC 5063 JCLJ 
3A 0557-383 @ 
IRAS 18325-5925 M 
IRAS 19580-1818 (42 j 
3C 324 @ 
3C 256 [44] 
3C 356 || 
F J084044.5+. . . 
F J155633.8+... 
3CR 68.1 fH 
QSO J2359-1241 
3C 234 M 



ll 



4C 40.36 [5j 
4C 48.48 @ 
IAU 0211-122 f0| 
IAU 0828+193 



L eff (Gpc) 



0.04 
0.12 
0.07 
0.14 
2.44 
3.04 
2.30 
2.49 
2.75 
2.48 
2.01 
0.61 
3.35 
3.40 
3.40 
3.53 



0.56 - 2.8 
2.2 - 8.5 
1.0 - 4.9 
1.8 - 9.3 
82 - 180 
110 - 220 
78 - 170 
88 - 170 
99 - 160 
84 - 180 
110 - 120 
55 - 81 
120 - 260 
120 - 260 
120 - 260 
130 - 270 



logio cr 



-30.8 
-31.2 
-31.0 
-31.0 
-32.3 
-32.4 
-32.3 
-32.4 
-32.4 
-32.4 
-31.2 
-31.7 
-32.4 
-32.4 
-32.4 
-32.4 



Table 2. Source data for polarization constraints. 

Table 2 lists 16 sources with published values of ip. The 
second column of the table displays the effective distance 
L e ff traveled by the light. The third column provides the 
range of wavelengths for which data are used. In fitting 
to 6ip, we choose to set ipo equal to the mean value of the 
measured ip. For each source, \q and Ao are fitted to the 
data. These angles can be thought of as the two degrees 
of freedom needed to describe the unknown polarization 
at the source. Ideally, at this point the data would be 
fitted to all the sources simultaneously. However, since 
k a has 10 elements and each source introduces two ad- 
ditional parameters, this would be involved. Instead, we 
examine each source individually and look for the desired 
wavelength dependence. 

Adopting the same analysis strategy as described in 
Ref. Jl^] yields the bounds for each source listed in the 
last column of Table 2, which can be combined to yield 
a bound of 



\k a \ < 2 x 10 



-32 



(27) 



in the Sun-centered celestial equatorial frame, at the 90% 
confidence level. 



IV. LABORATORY TESTS 

The Lorentz-violating electrodynamics predicts shifts 
in cavity- resonance frequencies, which offers the oppor- 
tunity for sensitive tests of Lorentz symmetry in labora- 
tories on the Earth and in space. This section presents a 
general framework for the analysis of such experiments. 
We begin with some general considerations and then sep- 
arately consider in more detail the cases of optical cavities 
and microwave cavities. 



A. General Considerations 

Many tests of special relativity search for variations 
in some observable that might arise from the rotation 
or boost of the apparatus due to the motion of the 
Earth. Lorentz-violating theories predict periodic vari- 
ations at multiples of the Earth's sidereal or orbital fre- 
quencies. For example, high-sensitivity measurements of 
coefficients in the fermion sector of the standard-model 
extension have been performed by comparing two clocks 
as the Earth rotates jy]. The clocks are typically the 
frequencies associated with specific Zeeman atomic tran- 
sitions, and the standard- model extension predicts vari- 
ations in these frequencies with the orientation of the 
apparatus and hence with the Earth's rotation. Similar 
tests could be performed in space, with the frequency 
variations depending on the orbital and rotational prop- 
erties of the spacecraft p0[ . 

Resonant cavities can also serve as clocks, and they can 
be used in clock-comparison experiments to test prop- 
erties of electromagnetic fields instead of atomic tran- 
sitions. In particular, clock-comparison experiments of 
this type can be used to probe the photon sector of 
the standard-model extension. One relevant issue in the 
analysis of these experiments is establishing the trans- 
formation between the laboratory frame and a standard 
celestial frame. Another is the determination of the pre- 
dicted frequency shifts. In this subsection, these issues 
are addressed in a general context. 



1. Generic laboratory experiment 

Consider a general laboratory-based experiment mea- 
suring some electrodynamic observable O. Typically, the 
constitutive relations (Q) change the observable from its 
conventional value Co- We consider a change SO, taken 
to be linear in the matrices kde, k-hb, and kqb = 



-K 



HE- 



In a frame fixed to the laboratory, SO can be 



written as 



SO = (M DE )£ h (K DE )£ h + (M 



\jk / \jk 
HB)ub\ K HB)ub 



where (Mde)^, (Mhb)^, and {M D B)ub are 
experiment-specific constant matrices determined by the 
apparatus. The symmetries of the n matrices can be im- 
posed on their M. counterparts when convenient. 

Due to the orbital and rotational motion of the Earth 
or the space platform, the laboratory cannot be con- 
sidered an inertial frame. As a result, the laboratory- 
frame coefficients (kde)^, (ksb)]^, and (kdb)^ vary 
in time. We can exploit the induced variation in O by 
searching for periodic fluctuations in O at the relevant 
frequencies. A measurement of this type of variation 
would be a signal for Lorentz violation. 

To determine the dependence of the periodic variation 
on the coefficients (kp)^^, we seek an expression sim- 
ilar to Eq. ( |2£| ) in an inertial frame. A suitable choice 
for a standard inertial frame is the Sun-centered celestial 
equatorial frame defined in Sec. Ill A. The coefficients 



,JK 



for Lorentz violation in this frame, (kde) > , (khb) 
and {kdb) JK , are constant. 

The observer Lorentz transformation between the two 
frames can be used to relate the corresponding two sets 
of k matrices. Since the velocity of the Earth with re- 
spect to the Sun is ~ 1CP 4 , it suffices for our purposes 
to construct the transformation to leading order. At this 
order, the Lorentz matrix implementing the trans- 
formation from the Sun-centered frame to the laboratory 
frame is 



A 



' T = l, A° 7 = -/? J , V T = -(R-py, A' ; //"'. 



(29) 



where f3 is the velocity of the laboratory with respect to 
the Sun-centered frame and R? J is the spatial rotation 
from the Sun-centered frame to the laboratory frame. 
Some calculation shows that the induced transformation 
between the k matrices is 

(A-D-EJlab — J ( K DE) -Jl (KDB) 

(KffBJlab - 1 \ K HB) — J-i (Kdb) , 

{KDB)u h =T^ {k DB ) J 



+T 1 kjJK (K DE ) JK 



rpjkJKf sJK 

■ J-l (KHB) , 



(30) 



where 



T; 



■jkJK 







W'R 



J nkK 



T 



jkJK 



R jP R kJ e KP Q/3 Q (31) 



(Mus)f ab (KDB)f ab , 



(28) 



The tensor Tq is a rotation, while T\ is a leading-order 
boost contribution. Although the contributions involv- 
ing Xi are suppressed by /?, they access distinct combi- 
nations of coefficients and can introduce different time 
dependence, which may lead to fundamentally different 
tests. 

To apply Eqs. (|30| ) and (^lj) , the laboratory frame must 
be specified. Appendix |C| defines our standard Earth- 
based and space-based frames and establishes the trans- 
formations from these to the reference Sun-centered ce- 
lestial equatorial frame. 



2. Cavity experiments 

Two classes of cavities are of interest in the present 
context: optical cavities, for which the wavelength of the 
light is much smaller than the cavity size, and microwave 
cavities, for which the wavelength and cavity size are 
comparable. In both cases, the interesting quantity is 
the fractional resonant- frequency shift 8v/v. 

For a given cavity, let Eq, Bo, Do, Ho be the fields 
associated with a conventional mode of resonant angular 
frequency uio- Nonzero kp coefficients can perturb these 
resonance fields. Let E, B, D, H be the perturbed fields 
for the resonant mode in the presence of Lorentz viola- 
tion, and let 8v — 8lj/2tt represent the change in the 
resonant frequency relative to the conventional case. A 
manipulation of the Lorentz- violating Maxwell equations 
then yields the fractional resonant-frequency shift as 



Su 
v 



d 3 x(E* -D + H*-B] 



< J d 3 x ^E* ■ D - D* ■ E - B* ■ H + H* ■ B 

-iw^V ■ (Ho x E - E* x H)j , (32) 

where the integrals are over the volume V of the cav- 
ity. This equation holds for any harmonic system, even 
for large differences between the conventional and per- 
turbed modes. Note that the divergence term results in 
a surface integral over the boundary of V. 

For the application to Lorentz violation, the perturbed 
modes are expected to differ only slightly from the un- 
perturbed ones. Also, the boundary conditions can rea- 
sonably be taken such that the divergence term in Eq. 
( |32l ) vanishes. The point is that, for leading-order ef- 
fects, we can approximate the cavity as lossless and ide- 
alize the surface of the cavity as a perfect conductor. The 
boundary condition of vanishing surface tangential elec- 
tric field Eo follows as usual from the Faraday equation 
V x £ + fl B = and the vanishing of Eo inside the 
conductor. The latter can be regarded as a consequence 
of the Lorentz force law. To determine the tangential 
perturbed field E on the cavity surface, we note that 
Lorentz violation in the photon sector leaves the force 
law unaffected. Disregarding for simplicity any effects 
on the force law arising from Lorentz violation in the 
fcrmion sector of the standard- model extension, which in 
any case would be expected to enhance a signal, it fol- 
lows that the tangential component of E also vanishes on 
the surface. With these boundary conditions, the normal 
component of (Hq x E — Eq x H) is zero at the surface 
of the cavity. 

For leading-order effects of Lorentz violation, it suffices 
to expand the remaining terms of Eq. ( ^2|) in the coef- 
ficients {kp) K \^. If the cavity is void of matter, then 



Do = Eo, Ho = Bo, and the constitutive relations (Q) 
yield the approximate relations 

D — E ~ kde ■ Eo + kdb ■ Bo, 

H - B ~ k he ■ E + k hb ■ Bq. (33) 

If the cavity contains matter, we adopt instead a general 
linear relation between the unperturbed fields {Do, Ho) 
and (Eo, Bq) and assume for simplicity a lossless medium. 
In either case, we find that the leading-order fractional 
frequency shift becomes 



1 



5v 



v A{U) 



j d 3 x (^Eo ■ kde ■ E - B*o ■ khb ■ Bo 



-2Re (E* ■ K DB -Bo)), (34) 



where (U) = J v d 3 x (E ■ D^ + B ■ #o)/4 is the time- 
averaged energy stored in the unperturbed cavity. Note 
that bvjv is real, reinforcing the argument that the vac- 
uum is lossless and indicating that the Q factor of the 
cavity remains unaffected by Lorentz violation at lead- 
ing order. 



B. Optical cavity experiments 

Among the classic tests of Lorentz invariance are the 
Michelson-Morley |l| and Kennedy-Thorndike Q exper- 
iments. Both concern the speed of light, with the former 
searching for spatial anisotropy and the latter seeking 
dependence on the laboratory velocity. The standard- 
model extension can be used as a general framework for 
analyzing these experiments. In this section, we consider 
modern versions of these tests that use optical cavities to 
achieve improved sensitivities 17- Tsf| . 



1. Theory 



We can use the results in Sec. [VA to obtain an ex- 



pression for the fractional frequency shift 8v/v arising 
from Lorentz-violating effects in an optical cavity. The 
idea is to regard the cavity as two parallel reflecting pla- 
nar surfaces with plane waves propagating between them 
normal to the surfaces, and then to apply Eq. (|34|). 

The resonant modes of optical cavities can be regarded 
as standing waves. For simplicity and definiteness, we 
suppose the unperturbed cavity contains a medium hav- 
ing transverse relative permittivity e and relative perme- 
ability (i = l, with the case e = 1 corresponding to a 
cavity void of matter. As usual, the unperturbed fields 
can be taken as 



E (x) = E cos{ljoN ■ x + (f)e- luJot , 
Bo{x) = iy/eN x E sin(w V / eA r ' 



-iuot 



(35) 



where N is a unit vector pointing along the length of the 
cavity, (ft is a phase, and Eq is a vector perpendicular to 
N that specifies the polarization. The conventional reso- 
nant frequencies are given by u>q — irm/y/el, where m is 
an integer and I is the length of the cavity. 

Substitution of Eq. ( p5|) into Eq. yields the desired 
result for the fractional frequency shift: 

Sv 1 

— (N x Eq) ■ {n HB ) l&h -(NxE )]. (36) 

This expression for the fractional frequency shift is also 
obtained in an alternative approach from a different phys- 
ical perspective, as described in Appendix 

The laboratory- frame matrices Mi^b introduced in Eq. 
(p8|) can be extracted from Eq. (|36|). We find 



k Re(EZy(E y 

( Mde ^ = 



2e\E Q \ 

Re (N x Eo*) j (N x E ) k 



2|£o 12 



(Mdb)^ = 



(37) 



These equations show that in the presence of Lorentz 
violation the frequency of an optical-cavity oscillator de- 
pends both on the orientation of the cavity and on the 
polarization of the light with respect to the laboratory 
frame. 

To analyze an experiment with an optical cavity, 
one can now proceed as follows. First, determine the 
laboratory- frame matrices Mi^b from the apparatus by 
applying Eq. (fTif). These matrices are constant if the 
cavity is fixed in the laboratory but vary with time if 
the cavity is rotated in the laboratory. Next, relate the 
laboratory-frame matrices Ki a b to those in the reference 
Sun-centered frame using the transformation ( |30| ) and 
the material in Appendix ^|. The time dependence of 
the cavity resonant frequency can then be calculated us- 
ing Eq. ( |36] ) or equivalently Eq. ( p8| ) . Finally, the ampli- 
tudes and phases of particular harmonics can be obtained 
and compared to the experimental data. 

As an illustration of the analysis procedure, consider 
laser light incident on a cavity positioned horizontally in 
an Earth-based laboratory, with the light linearly po- 
larized along the z axis. Denote by 9 the angle be- 
tween the x axis and the cavity orientation. Then, 
N = (cos 6, sin 6,0), and in the laboratory frame the frac- 
tional frequency shift becomes 

Sv 

V 

-h(KHB)l* h shi29 

- ("*u)2b] cos2tf. (38) 



The next step is to transform this result to the Sun- 
centered celestial equatorial frame. Using Eqs. ([50|) and 
(pl|), the fractional frequency shift takes the form 



— = A 
v 



B sin 29 + C cos 29, 



(39) 



where 



A = A + A\ sinw®T ffi + A 2 cos uj®T s 
+A 3 sin 2w e T e + A 4 cos 2w®T e , 

B = B + Bi sinu; e r e + B 2 cosw e T(j 
+B 3 sin 2w e T® + B4 cos 2w©T® , 

C = Co + C\ sinw e T© + C2 coswqT® 
+C3 sin 2w e T e + C 4 cos 2w e T e . 



(40) 



The quantities ^0,1,2, 3, 4, -80,1,2, 3, 4, and 60,1,2,3,4 are lin- 
ear in the coefficients for Lorentz violation and depend 
on the colatitude x- They are given explicitly to order 
f3 in Appendix |]. Note that the coefficient k tv appears 
only in Aq, resulting in a constant frequency shift. It fol- 
lows that sensitivity to k tI is suppressed by at least two 
powers of P in this experiment. 

The analysis could now proceed along several lines. 
One possibility is to adopt the birefringent constraints 
(p7|). The expressions in Appendix Je] can then be sim- 
plified by setting (k e+ ) JK = (k -) = 0. This shows 
that the eight coefficients directly accessible 

through fitting the measured frequency shift to Eq. (|4"c|). 
Another possibility is to include all coefficients in the 
analysis. This would provide a direct laboratory check on 
the birefringence results. Although in practice the sen- 
sitivity is much reduced, the systematics of laboratory- 
based experiments are fundamentally different from those 
in cosmological tests and so this check is worthwhile. 
We remark that the isolation of specific coefficients for 
Lorentz violation can be aided by considering different 
experimental configurations. These include adopting a 
different polarization and rotating the apparatus in the 
laboratory, which produces a time dependence in 9. 



2. Experiment 



A modern Michelson-Morley experiment with sensitiv- 
ity to a fractional frequency shift 8v/v of about 10~ 13 
was performed by Brillet and Hall [ pT[. A similar sensi- 
tivity was achieved by Hils and Hall |18|] in a Kennedy- 
Thorndike experiment, recently repeated using a cryo- 
genically cooled cavity by Braxmaier et al. |l9) . These 
experiments compare the fractional frequency shifts be- 
tween two lasers. One laser is stabilized to a molecular 
transition and serves as a reference frequency. A portion 
of the light from the second laser is sent into one end of an 
optical cavity, and the light emerging at the other end is 
used to tune this laser to the cavity resonant frequency. 



The remaining light from the second laser is combined 
with the light from the reference laser, and the beat fre- 
quency is measured. In the classic analysis, the frequency 
of the reference laser is independent of violations of spe- 
cial relativity while the frequency of the cavity-stabilized 
laser depends on the speed of light along the length of 
the cavity. 

The Brillet-Hall experiment studies spatial isotropy by 
seeking changes in the beat frequency as the cavity is 
rotated in the laboratory with a period of about 10 s. 
The vector Fourier amplitude is measured at twice the 
cavity rotation frequency. In the present context, if we 
suppose for definiteness a vertical laser polarization as 
in the previous subsection, this experiment offers sen- 
sitivity to the quantities B and C in Eq. (|9|) through 
the time dependence of 9. The reported fractional fre- 
quency shift is 1.5 ± 2.5 x 10~ 15 . The analysis yielding 
this bound supposes a signal at 2w® and averages over 
several days of data. Within the framework leading to 



Eq. (40), this bound would translate to a constraint on 
a particular combination of the coefficients -B3, B4, C3, 
C4 at the level of about a part in 10 15 . 

A complete dataset of the type taken in this experi- 
ment could be analysed using Eq. ( p0| ) to extract several 
different measurements of combinations of B n and C n . 
For example, consider the one-day dataset displayed in 
Fig. 2 of Ref. In this dataset, no variation is seen in 
the frequency above the level of \J B 2 + C 2 < 4 x 10~ 13 . 
In fact, these data exhibit a one-day signal involving a 
roughly constant Fourier amplitude of about 2 x 10~ 13 
with nearly constant phase, attributed to a slight tilt 
in the rotation platform. Since a nonzero value of 
y Bq + Cq would produce a similar signal, compelling 
measurements of Bq, Co via this method appear prob- 
lematic. However, bounds on combinations of the quan- 
titites B n , C n for n ^ could be extracted by studying 
the behavior of the data at both the sidereal frequency 
and its harmonic 2lo@. As can be seen from the ex- 
pressions in Appendix the quantities B n , C n involve 
unsuppressed combinations of the coefficients k e+ , k e _ 
for Lorentz violation, along with combinations of the co- 
efficients k 0+ , suppressed by one power of the veloc- 
ity. It therefore appears feasible to perform a systematic 
analysis of a complete dataset in a Michelson-Morley ex- 
periment with an optical cavity to measure combinations 
of the coefficients k e+ , K e _ with a sensitivity of order 
2Q-i4±i an( j combinations of k +, k - with a sensitivity 
of order 10 _10±1 . 

The Hils-Hall experiment seeks changes in the beat fre- 
quency as the velocity of the laboratory varies with the 
Earth's rotation. The analysis assumes that experiments 
of the Michelson-Morley type exclude observable sensi- 
tivity to the orientation 9 of the cavity. In the context of 
Eq. (^), this corresponds to assuming negligible B and 
C terms. The Fourier amplitude at the sidereal frequency 
is analyzed, obtaining a bound of 2 x 10~ 13 at the 90% 



confidence level. With the configuration leading to Eq. 
((40|), this bound constrains the combination \J A\ + A 2 . 

Since at present many combinations of B and C remain 
unconstrained, the assumption of negligible B, C terms 
is undesirable in the analysis of Kennedy-Thorndike ex- 
periments. If this assumption is relaxed, the Fourier am- 
plitude at the sidereal frequency contains contributions 
from A\ , A2 , B\ , B2 , C\ , C2 ■ It should therefore be pos- 
sible to measure combinations of the coefficients k e+ , k e - 
at the level of about 10 -13 and combinations of the coef- 
ficients k Q+ , k ~ at the level of about 10 -9 . A complete 
analysis could also analyze the second Fourier amplitude, 
which would provide a measurement of a combination of 
^3, A4, B 3 , B 4 , C3, C4. 

The analysis of Braxmaier et al. focuses on a varia- 
tion in 8v/v with the orbital motion of the Earth. The 
assumption of negligible £?, C terms is again made. The 
data are averaged daily, and a bound on the fractional 
frequency shift of 4.8 ± 5.3 x 10~ 12 is obtained. In the 
context of Eq. (fl0"|), the analysis restricts the sensitiv- 
ity to Aq. Using Eq. (El) in Appendix || it then follows 
that the reported constraint corresponds to a bound on a 
combination of k a+ , k a ^ at the level of about 10 -8 . Sen- 
sitivity to Co could also be obtained if B and C terms 
were included. Note that the polarization chosen in de- 
riving Eq. ( [40| ) implies the coefficient Bq is independent 
of /?9 and hence cannot be extracted. 

The above discussion shows that many interesting 
prospects remain for measurements of the coefficients k 
using optical cavities. Note that the published exper- 
imental analyses to date are each sensitive to different 
combinations of coefficients for Lorentz violation. A sys- 
tematic analysis could in principle provide sensitivity to 
all the coefficients k e+ and K e _ at the level of about 10~ 13 
or better and suppressed sensitivity to the coefficients 
k a+ and k d _ at the level of about 10~ 9 or better. Note 
also that an analysis along the above lines could readily 
be applied to space-based tests involving optical cavities, 
including ones on the ISS or in dedicated missions such 
as the proposed OPTIS experiment |^l| . Some potential 
advantages of space-based tests are discussed below in 
the context of experiments using microwave cavities. 



C. Microwave cavities 

Microwave-cavity oscillators are among the most sta- 
ble clocks, and as such they offer interesting prospects 
for Lorentz tests. In particular, there has recently been 
renewed interest in superconducting cavity-stabilized os- 
cillators as clocks for use on the ISS [52] |. Superconduct- 
ing cavities made of niobium have achieved Q factors of 
10 11 or better, and frequency stabilities of 3 x 10 -16 have 
been demonstrated. In this section, we focus on pertur- 
bations of microwave-cavity resonant frequencies arising 
from the coefficients {kp) K \^ v for Lorentz violation. 



1. Theory 

Equation (|34|) can be applied to obtain the fractional 
resonant-frequency shift 5v/v for a superconducting mi- 
crowave cavity of any geometry. The highest Q factors 
have been demonstrated in cylindrical cavities with cir- 
cular cross section, so we focus on this case. For simplic- 
ity, we suppose the cavity contains a medium of relative 
transverse permittivity e, relative axial permittivity e', 
and relative permeability fi = 1. The vacuum case is 
recovered as the limit e = e' = 1. 

The invariance of the cylindrical geometry under a par- 
ity transformation suggests the matrix Mdb vanishes, 
since the coefficients {kdb)^ ar e parity odd. Also, the 
cavity is invariant under rotations about the symmetry 
axis, so in a cavity frame with 3 axis along the symmetry 
axis we expect the rotational symmetry to imply diagonal 
matrices (Mde)cav and (Mhb)c&v with equal {11} and 
{22} components. Indeed, for both TE m „ p and TM mnp 
modes, we obtain 



Si; 



(Mde)H v [( k de)H v + (kde)H v 
+ (M D e) 3 c I v {kde) 3 c 



33 
cav 
11 

cav 



-(Mhb)H v (khb)H v 



22 i 
cav J 



(41) 



in the cavity frame. 

For the TM mnp modes, some calculation reveals that 
the nonzero elements of the M matrices are: 



(M DE )ll = (M 
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e'{TipRf 
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cav 
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2 e' 2 (7rpR) 2 + ee'(x mn d) 2 - 



(M--^ 22 



(42) 



where R and d are the radius and length of the cavity, 
and where x mn is the nth zero of the mth-order Bessel 
function J m {x). The corresponding results for the TE mnp 
modes are 



(M D e)H v - (M D e)H v - -j-, 



(M 



HB 



(Mhb) 
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(7rpR) 2 + (x' mn d) 2 ' 



(Mhb)11 



(npR) 2 + (x' mn d) 2 ' 



(43) 



where x' mn is the nth zero of the derivative of J m {x). 
Note that taking the optical-cavity limit p — > oo of any 
TM mnp or TE mnp mode yields a result identical to that 
obtained by averaging over all optical-cavity polariza- 
tions in Eq. @. 

For practical applications, it is useful to generalize Eq. 
([il]) to the case where the cavity is arbitrarily oriented 



in one of the standard laboratory frames introduced in 
Appendix |^. In the laboratory frame, denote the com- 
ponents of a unit vector parallel to the symmetry axis of 
the cavity by . The fractional frequency shift is then 
found to be 



v 



= {M D E)H v (KDE)il h + {M H B)ll v (KHB) J d h 



-[(M DE )f av - (M DE )ll v ]NJN k (K DE )i* h 
i(M HB ) 3 c l - (M H B)ll v ]NiN k ( KHB )Q h . 



(44) 



This implies the relationships 
(M DE t h = [(Mdb) 



{M HB ){* b = [(Mhb)Z - (M HB ) cav 
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\33 



{M DE )HWN k 



]N 3 N k 



-(Mhb) 1 ^ 



0. 



(45) 



Using Eq. (44) and the transformation (p0|), we can 
write the fractional frequency shift for a general mode 
in terms of coefficients for Lorentz violation in the Sun- 
centered celestial equatorial frame. To order /3, we find 



5v 



-\N J N k R jJ R kK ^- 



(H e ,) jK -\(M 2 + 3M 3 )H tl . 

-\{M^ k /M 2 + N 1 N k )W J R kK e JPQ [3 Q (k ,) KP . 

(46) 

In this equation, we define the quantities 
Mi = -A[{M DE )\ 
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(47) 



which depend on the cavity mode and control the linear 
combinations 



(h e ,) JK = M 1 {k e+ ) JK + M2{k e -) JK , 

(k ') JK = Mi(k -) jk + M 2 (k 0+ ) JK 



(48) 



of coefficients for Lorentz violation. These equations re- 
veal that the sensitivity of experiments with microwave 
cavities to Lorentz violation varies with the mode and 
with the permittivity of the medium in the cavity. As 
before, to this order htr contributes only to an unobserv- 
able constant frequency shift. 

As an illustration, consider a cavity void of matter and 
operated on the fundamental TMoio mode, as planned 
for some space- and ground-based experiments. For this 
case, we find Mi — 3, M 2 — M3 — 1, and the fractional 
frequency shift (46) becomes 
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-\N3N k [2{K DE )£ h + {K HB )£ h 

-8 3k {nHB)u h ] 
\N3N k W J R kK [Z{k e+ ) JK + (£ e _) JA '] 

.1 + N^N k )W J R kK e JPQ f3 Q 

x[3(£ _)*" p + (~ Ko+ ) KP ]~k tr . (49) 



The observable shift depends on the traceless symmetric 
matrix combination 3(k e+ ) JK + (k e ^) JK and the trace- 
less matrix combination 3(k ^) JK + (k +) jk . The first 
of these contains five linearly independent combinations 
of the 11 parity-even coefficients for Lorentz violation, 
while the second contains all eight parity-odd coefficients. 
Note that certain harmonics may be sensitive to smaller 
subset of these 13 quantities. For example, for a fixed 
Earth-based cavity, if Cj represents the rotational axis of 
the Earth's revolution, then the sidereal harmonics are 
insensitive to the component of 3(k e+ ) JK + (k e _) JK pro- 
portional to Cj j Cj k , reducing the number of combinations 
to 12. 

If an Earth-based experiment is performed over a pe- 
riod of time AT exp short compared to the orbital pe- 
riod of the Earth, then the velocity (3 is roughly constant 
and the experiment is sensitive primarily to the four lin- 
ear combinations corresponding to the vector amplitudes 
of the two harmonics. To acquire sensitivity to other 
combinations, the Earth-based experiment could either 
be repeated several times during the year, or the cav- 
ity could be rotated in the laboratory. In contrast, for 
a satellite-based experiment, perturbations cause the or- 
bital plane and hence the analogue of Cj to precess with 
time. Also, the smaller orbital period implies different 
harmonics and access to more coefficients for the same 
AT exp . As a result, if the experiment is performed two or 
more times with significantly different orbital planes, all 
13 combinations of coefficients can be accessed through 
the orbital- frequency harmonics. 

Equation (Eh can be adapted to either a space-based 
or Earth-based experiment and, if necessary, to the case 
of a rotating cavity. In the remainder of this section, 
we offer some remarks about possible experiments with 
microwave cavities on the ISS and on the Earth. 



2. Space-based experiment 

The construction of the ISS offers the possibility of 
performing Lorentz tests in low Earth orbit. Of par- 
ticular relevance in the present context is the SUMO 
experiment [52] , which plans to use superconducting 
microwave-cavity oscillators as clocks on upcoming ISS 
flight missions. 

The ISS operates in several different flight modes, 
which correspond to different laboratory configurations 
in the Sun-centered celestial equatorial frame. Each flight 



mode therefore involves different transformations (|3C 
which could lead to different sensitivities to the Lorentz- 
violating coefficients. If, for example, the ISS orientation 
were fixed in Sun-centered frame, corresponding to no ro- 
tation during an orbit, then the signal would involve only 
boost-dependent terms with a 92-minute period. For def- 
initeness and simplicity, we focus here on a flight mode 
with the ISS z-axis aligned along its orbital velocity with 
respect to the Earth. This corresponds to the standard 
laboratory frame introduced in Appendix 

For a microwave cavity with fixed orientation N in this 
ISS laboratory frame, several harmonics could be stud- 
ied. The fractional frequency shift 5v/v varies with the 
Earth's orbital frequency, the ISS orbital frequency uj s , 
and the ISS orbital-precession frequency. The most in- 
teresting of these is likely to be the highest frequency, 

LU S . 

In practice, the fractional frequency shift may be mea- 
sured relative to another oscillator clock via the beat fre- 
quency of the combined signal. The reference clock could 
be a different physical system, such as a hydrogen maser 
or atomic clock, which could conveniently be operated 
on a transition known to be insensitive to Lorentz viola- 
tion ptJ. A comparison of two microwave cavities could 
also be used. For example, SUMO may involve a pair 
of cavities oriented at right angles to each other. The 
observed signal would then depend strongly on the ori- 
entation of the pair in the ISS frame. Thus, at leading 
order in /3, a cavity oriented with N perpendicular to the 
orbital plane is insensitive to the parity-even coefficients 
for Lorentz violation, and only /^-suppressed parity-odd 
terms appear in the frequency shift. In contrast, a cav- 
ity positioned with N in the orbital plane maximizes the 
sensitivity in the second harmonics of lo s , while one with 
N at 45° from the orbital plane maximizes the first har- 
monics. These results can be obtained directly from Eq. 

For a pair of identical cavities, the variation in the beat 
frequency takes the general form 

Vbeat _ SVI _ Sf 2 
V V V 

= A s sin U! S T S + A c cos uj s T s 

+B S sin 2w s T s + B c cos 2lu s T s + C, 



(50) 



where A s , A c , B s , and B c are four linear combina- 
tions of the coefficients {kp) K \^ l , for Lorentz violation. 
These combinations depend on the orientations Ni, N 2 
of the cavity pair and on the orientation of the orbital 
plane with respect to the Sun-centered celestial equato- 
rial frame. The precession of the ISS orbit slowly changes 
the four combinations, allowing access to more coeffi- 
cients. Typically, the combinations are rather cumber- 
some. Appendix [f] contains their explicit form for a 
maximal-sensitivity case, for which N\ = (0,0,1) and 
N 2 = (1,1,0)/V2. The expressions involve the linear 



combinations (fl8|), which hold for an arbitrary mode and 
arbitrary permittivities e, e'. 



3. Earth-based experiment 

For an Earth-based experiment with a cavity pair fixed 
in the laboratory, the dominant frequency is the Earth's 
sidereal frequency w©. The equivalent of the 1SS orbital 
plane in this case is the plane in which the laboratory 
moves, which parallels the equatorial plane at the lati- 
tude of the laboratory. As before, the configuration of 
maximum sensitivity has one cavity in this plane and 
the other at 45° to it. However, for laboratories located 
in middle latitudes, it suffices to orient one cavity hor- 
izontally in the east-west direction and the other either 
vertically or horizontally in the north-south direction. 
The east-west cavity is then maximally sensitive to the 
second harmonics, while the north-south cavity is near 
maximal sensitivity to the first harmonics. The latter 
are proportional to cos 2%, so for colatitudes in the range 
30° < x < 60° there is at most a 14% reduction in sensi- 
tivity. 

For definiteness, we consider the configuration with 
the second cavity oriented vertically in the laboratory. 
The laboratory-frame orientation vectors are then N\ = 
(0,1,0) and N 2 = (0,0,1). Paralleling the discussion 
leading to Eq. (|5(]) , we write the fractional beat frequency 
due to Lorentz violation as 

= „4© s sinu;©T© + .4© c cosw e T e 

v 

+B® S sin 2w©T© + B @c cos 2w e T ffi + C©. (51) 
At first order in /?, we find 
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where the convenient combinations (48) have been 
adopted. As before, these equations are valid for any 
specific mode and for arbitrary permittivities e, e'. 

A Lorentz-violating signal would thus manifest itself 
as a sidereal variation in the fractional beat frequency 
according to Eq. (|5l|). At zeroth order in f3, the cor- 
responding amplitude associated with this variation is 
constant and determined by the four parity-even coeffi- 
cients (k e ,) xz , (k d ) YZ , (k e ,) XY , and (R e >) xx -{R e >) YY ■ 
These linearly independent combinations of (kF) K \^i> re- 
main unmeasured to date. When the first-order terms 
in (3 are included, the amplitudes also contain harmonics 
at the Earth's orbital frequency fi©. The resulting vari- 
ations depend on the eight parity-odd coefficients k >. 
The three of these represented by k 0+ have yet to be 
measured. 

The above experiment provides access to the specified 
parity-even coefficients at the level of the cavity stabil- 
ity, which for available microwave cavities could be at the 
order of 10 -13 or better. The (3 suppression reduces the 
sensitivity to the parity-odd coefficients to about 10~ 9 or 
better. The experiments can be performed on any cavity 
mode and for cavities with or without matter. For ex- 
ample, as can be seen from the explicit expressions (BH), 
a pair of sapphire-filled cavities (e ~ 9.5, e' ~ 11.5) oper- 
ated on a whispering-gallery mode J53| offers sensitivity 
to linear combinations of coefficients for Lorentz viola- 
tion that differ from those of a pair of vacuum cavities 
operated on the fundamental TM010 mode. In fact, for 
any specified combination of cavities with known fields, 
the matrices (Mde) C !xv, {Mhb)c&-v in Eq. (41) can be 
determined and hence the fractional beat frequency can 
be obtained as above. Note also that other coefficients 
could be accessed by placing the cavity pair on a rotat- 
ing turntable, which would also allow a dataset to be 
obtained in weeks or days rather than months. 



D. Cavity Deformation 



In this remaining subsection, we offer a few remarks 
concerning the possibility that Lorentz violation might 
alter atomic binding forces and hence the structural prop- 
erties of matter. In particular, one resulting effect rele- 
vant to cavity experiments might be a deformation of 
the cavity, which could change the resonant frequency. 



The issue is whether possible deformations arising from 
Lorentz- violating effects in the standard-model extension 
could cancel the predicted effects. The differences be- 
tween the transformation properties and nature of the 
various coefficients for Lorentz violation make it un- 
likely that coefficients other than (fcj?) re A J ui/ could cause 
complete cancellation of the signals discussed above, so 
a more interesting question is whether the coefficients 
(kFjnXtif could alter the dimensions of a cavity so as to 
offset completely the predicted signals. 

Any leading-order modifications to atomic and molec- 
ular binding forces arising from nonzero (kF) K \^v are ex- 
pected to come from modifications to the Coulomb po- 
tential. The form of the Gauss law (||) in the presence of 
Lorentz violation due to (/cf) k aa"/ implies that the mod- 
ified Coulomb potential for a point charge q is 
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The leading-order effects on the physical dimensions of 
the cavity are therefore expected to depend only on the 
matrix kjje- We could account phenomenologically for 



such effects by adding a term to Eq. ( |28[ ) of the form 

jk 



M 



4jk 



lab' w here the constant matrix M matteI 



determined by the properties of the material from which 
the cavity is made. For example, in a simple ionic lattice 
model, M^ attcr depends on the charge of the ions, the 
lattice configuration, and the orientation of the cavity 
with respect to the lattice. 

For the optical and microwave cavities considered here, 
this extra term cannot completely cancel the predicted 
fractional frequency shifts. Although a partial cancella- 
tion might be possible in principle, it requires that the 



matrix M. 



oh 



takes a special form that is improbable 



in light of the complexity of the binding forces of solids. 



V. SUMMARY 

In this work, we studied the Lorentz-violating electro- 
dynamics derived from the renormalizable sector of the 
full Lorentz-violating standard-model extension. Some 



basic material is presented in Sec. II A, followed by a 
useful analogy and some definitions in Sec. II B. Sec. 



[I C discusses the connection to some test models. 

The bulk of the paper is devoted to tests of the Lorentz- 
violating electrodynamics and methods to measure the 19 
independent coefficients (fop^A/ii/ for Lorentz violation. 
We first consider astrophysical tests based on the predic- 
tion that the vacuum is birefringent. Theoretical issues 
perta ining to vacuum birefringence are discussed in Sec. 
|lll A| . One potentially observable effect is the dispersion 
of pulses over astrophysical distances. The constraint 
( pi| ) on (fcir)«AjL«/ from this effect is obtained in Sec. 



IIIB 



Another potentially observable effect arises in the com- 
parative spectropolarimetry of cosmological sources. The 



tight bound ( p7| ) on 10 of the 19 coefficients (A;f)kA/w is 
obtained in Sec. Ill C . 



The possibility of sensitive laboratory tests of Lorentz 



invariance is examined in Sec. IV A. A general framework 
for the analysis of both Earth-based and space-based ex- 
periments is provided. The analysis is applied to two 
types of cavity-stabilized oscillator experiments. In Sec. 
[VB , we consider optical-cavity experiments. High sensi- 
tivity microwave-cavity experiments are discussed in Sec. 
[V C| . We find that appropriate laboratory tests can ac- 
cess all 19 coefficients {kp) K \^ v . 

Table 3 summarizes the existing constraints. The 19 
coefficients {kp) K \^ v are represented by the matrices H e+ , 
K e _, k Q+ , k _, K tr defined in Eq. ([?]). The number of in- 
dependent components in each matrix is shown in the sec- 
ond column. The order of magnitude of the astrophysical 
bounds is shown in the third and fourth column. These 
bounds tightly constrain the 10 coefficients {kp) K \ fJil , con- 
tained in K e+ and ft G _. However, as indicated in the ta- 
ble by the notation n/a, the remaining coefficients are 
unobservable in astrophysical tests. In contrast, labo- 
ratory experiments with optical and microwave cavities 
can in princ iple access all the coefficients. As discussed 
in Sec. IV B| , several recent experiments with optical cav- 
ities |fl7Hl9[ offer sensitivity to a few of the coefficients 
at levels lying between about 10~ 8 and 10 -15 , but no 
definitive analysis has been performed. The matrices for 
which a few components are probably constrained in this 
way are indicated by the symbol * in the table. To date, 
no measurements of Lorentz violation using microwave 
cavities have been reported. 





Astrophysical Tests 


Laboratory Tests 


Coeff. 


No. 


Velocity 


Polarization 


Optical 


Microwave 




5 


-16 


-32 






K e — 


5 


n/a 


n/a 








3 


n/a 


n/a 






K — 


5 


-16 


-32 








1 


n/a 


n/a 







Table 3. Existing constraints. 



In conclusion, astrophysical observations place bounds 
on Lorentz violation in electrodynamics that are compet- 
itive with ones in the fermion sectors obtained by other 
means. Laboratory experiments are needed to comple- 
ment these measurements by spanning the allowed pa- 
rameter space in the photon sector, and the technology 
presently exists to perform them. These experiments of- 
fer a promising avenue to search for new physics lying 
beyond the standard model. 
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APPENDIX A: BIREFRINGENCE VECTORS 



For a distant source viewed from the Earth at dec- 
lination d and right ascension r, the direction of prop- 
agation towards the Earth can be written as p^ = 
(1 ; — cos d cos r, — cos d sin r, — sin d) . The matrix 



W K = 



sin d cos r 

sinr 
- cos d cos r 



sin d sin r 

— COST 

- cos d sin r 




(Al) 



implements the rotation between the primed frame and 
the standard Sun-centered frames. With this definition, 
the primed-frame basis vector e' 3 points from the source 
towards the Earth. The vectors e[ and e' 2 point south 
and west, respectively. 

Writing a sin £ and <rcos£ in terms of coefficients in 
the Sun-centered celestial equatorial frame gives 



a sinC - \{R 1J R 2K + R 2J R 1K )k J / Ku p„p v , 
acos£ = \{R lJ R lK - R 2J R 2K )k J / Kl/ p^. 



(A2) 



Note that £ is not a rotational scalar, unlike p and a. 



The rotation (Al) can be substituted in this result to 
yield cr sin £ and <rcos£ in terms of (fcir) K X/ 1 w in Sun- 
centered celestial equatorial coordinates. The relevant 
combinations of the (&f )«A/w are the 10 coefficients k a 
given in Eq. (|J) . It is convenient to express a sin £ 
and crcos£ as the scalar product of k a with two 10- 
dimensional vectors. Defining 

/ cos 2 d + cos 2 r — sin 2 d sin 2 r \ 
sin 2 d cos 2 r — cos 2 d — sin 2 r 
— 2 sin d sin r cos r 
— sin d sin r cos r 



sin d{sm 2 



— cos d sinr 
cos d cos r 

- sin d cos d cos r 
cos 2 d sin r cos r 

- sin d cos <i sin r 



(A3) 



we find 



— 2 sin <i sinr cos r 
—2 sin d sin r cos r 
~(1 + sin 2 d)(sin 2 r - cos 2 r) 
i (sin d + sin 2 r — sin 2 d cos 2 r) 
(1 + sin 2 d) sin r cos r 

— sin d cos <i cos r 

— sin <i cos d sin r 

cos d sin r 
sind(sin 2 r — cos 2 r) 
— cos d cos r 

asin£ = ^k a , 
crcosf = <; c a A: Q . 



(A4) 



(A5) 



APPENDIX B: POLARIZATION REVIEW 

Conventionally, polarization is defined by the behavior 
of the electric field vector . The polarization of a gen- 
eral plane wave can be described by an ellipse residing 
in the plane perpendicular to the direction of propaga- 
tion. In terms of the primed-frame variables introduced 
in Sec. Ill A and to leading order in (kp) K \^, this plane 



is spanned by the basis vectors e\ and e! 2 ■ The orientation 
and shape of the ellipse can be described by two angles, 
if) and x- The angle ip determines the orientation of the 
ellipse and is defined as the angle between the major axis 
of the ellipse and e[ . The angle x describes the shape of 
the ellipse and the helicity of the wave, and it is defined 
by x — i arctanr, where r is the ratio of the minor to 
major axes of the ellipse. 

In polarimetry, the ellipse is commonly parametrized 
using Stokes parameters. We define a Stokes vector by 



( S °,s)^(\E[\ 2 + 

2ReE[*E' 2l 2ImE[*E' 2 ) 
= s (1, cos 2% cos 2-0, cos2xsin2?/', sin2x). 



'2W \E[\ 2 ~\E' 2 \ 2 , 



(Bl) 



In the context of the discussion in Sec. 1IIC, the loss 



lessness of the vacuum implies that the Stokes parame- 
ter s° is unaffected at leading order by a relative-phase 
change. We therefore normalize to s° = 1 throughout. 
With s° = 1, each Stokes vector s represents a unique 
point on a two-dimensional sphere of unit radius, called 
the Poincare sphere. As illustrated in Fig. ^, 2ip and 2x 
are the angles that specify the position of s on this sphere. 
An arbitrary polarization is represented by a single point 
on the sphere. The points in the s 1 -s 2 plane represent all 
linear polarizations. The points in the upper hemisphere 
all represent elliptical polarizations of positive helicity, 
with the pole being the special case of circular polariza- 
tion. Similarly, the lower hemisphere represents polar- 
izations of negative helicity. 



circularly polarized 




. circularly polarized 
FIG. 2. The Poincare sphere. 



In Sec. |III Cj the quantity of interest is the polariza- 
tion change induced by the phase shift A(j> in Eq. (23). 
The effect of A</> on the Stokes vector s can be visualized 
in terms of motion on the Poincare sphere. Consider an 
arbitrary orthonormal elliptical basis {ei,£2}- The asso- 
ciated Stokes vectors sg t , sg 2 determine opposite points 
on the Poincare sphere. Decomposing a general electric 
field in this basis gives polarization components E n e~ l ^ n , 
n = 1,2, where E n and <p n are real. Examining the 
Stokes vector for this configuration shows that a change 
A (0i — 4>2) in the relative phase results in a right-handed 
rotation of the Stokes vector by the angle A(0i — 4>2) 
about the axis given by sg 1 — —sg 2 . 



APPENDIX C: STANDARD FRAMES 



cos x cosw®T® cosxsinw^Xg 
R jJ = I -sincj e T e cosw©T e 

sinxcosw e T e sinxsin^eTg; 




(CI) 



In this equation, j = x,y,z — 1,2,3 denotes an index 
in the laboratory frame, while J = X, Y, Z denotes an 
index in the Sun-centered frame. The Earth's sidereal an- 
gular frequency is a% ~ 27r/(23 h 56 min.), and x is the 
colatitude of the laboratory. The time T® is measured 
in the Sun-centered frame from one of the times when 
the y- and Y-axes coincide, to be chosen conveniently for 
each experiment. The time Tq therefore differs from the 
celestial equatorial time T by a constant shift for each 
experiment. 

The velocity 3-vector of the laboratory in the Sun- 
centered frame is 



sinllgT 
cos?7COsf2 ffi T 



cosw e T e 




(C2) 



Here, f2® and /3® are, respectively, the angular frequency 
and speed of the Earth's orbital motion. The quantity 
r) ~ 23.4° is the angle between the XY celestial equa- 
torial plane and the Earth's orbital plane. The speed 
@L = fffi^esinx < 1.5 x 10~ 6 is that of the laboratory 
due to the rotation of the Earth. 

The reader is warned that the standard laboratory 
frame defined above may differ from a frame fixed to the 
apparatus in the laboratory. For example, the apparatus 
rotates in the laboratory in some experiments considered 
here. Where confusion could occur, we distinguish with 
labels the quantities defined in the standard laboratory 
frame from those in the apparatus frame. 



This appendix defines our standard frames for Earth- 
and space-based laboratories. We provide the rotations 
and velocities used in transforming to the reference Sun- 
centered celestial equatorial frame, which is defined in 
section III A. 



1. Earth-based laboratory 



For a laboratory fixed to the surface of the Earth in 
the northern hemisphere, we choose the standard frame 
to have coordinates (t, x, y, z) such that the x-axis points 
south, the y-axis points east, and the z-axis points verti- 
cally upwards. With the reasonable approximation that 
the orbit of the Earth is circular, the rotation from the 
Sun-centered celestial equatorial frame to the standard 
laboratory frame is given by 



2. Space-based laboratory 

For our standard laboratory fixed to an Earth-orbiting 
space platform such as the ISS, we choose the z axis to be 
aligned with the velocity $ s of the satellite with respect 
to the Earth. The x axis is chosen to point towards the 
Earth. The y axis completes a right-handed coordinate 
system, thus directed along the satellite orbital angular 
momentum with respect to the Earth. 

The components of the matrix R? J describing the rota- 
tion from the reference Sun-centered frame to this stan- 
dard satellite frame are 



R 



IX 



cos a cos lu s T s + sin a cos £ sino; s T s , 
R 1Y = — sinacosuj s T s — cos a cos £ sinw s T s , 
R 1Z = -sinCsinw s T s , 
R 2X = sin a sin (, 



R 2Y = — cos a sin 
R 2Z = cosC, 

R? x — — cosasinw s T s — sinacosCcosw s T s , 
R 3Y = — sin a sin to s T s + cos a cos £ cos uj s T s , 
R 3Z = sinCcosw s T s . (C3) 

Here, j = 1,2,3 denotes an index in the satellite frame. 
The satellite orbital angular frequency is denoted u> s . 
The time T s is measured in the Sun-centered frame from 
a conveniently chosen time when the satellite crosses 
the equatorial plane, so the times T s and T differ by 
a constant for each experiment. Also, £ is the angle be- 
tween the satellite orbital plane and the Earth's equa- 
torial plane. For example, £ ~ 52° for the ISS. The 
quantity a is the azimuthal angle at which the orbital 
plane intersects the Earth's equatorial plane. The satel- 
lite intersects the equatorial plane twice per orbit, and 
a can be regarded as the angle between the X direction 
and a vector from the Earth's center to the point where 
the intersection occurs on an ascending trajectory. 

The velocity of the satellite with respect to the Sun- 
centered celestial equatorial frame is 

P — \ — cos77Cosn e T 
- sin r\ cos 

cos a sin ui s T s — sin a cos £ cos ui s T s 
sinasino; s T s + cos a cos Q cos uu s T s . (C4) 
sm. cos uj s T s 

The quantities are defined as before. The quan- 

tity p s is the speed of the satellite with respect to the 
Earth. For example, ft~3x 1CT 5 for the ISS. 



APPENDIX D: OPTICAL FREQUENCY SHIFT 

This appendix provides an alternative method to ob- 
tain the fractional resonant- frequency shift bvjv for op- 
tical cavities, given in Eq. (36). As before, the cavity is 
regarded as two reflecting parallel planar surfaces sepa- 
rated by a distance I. For convenience, we let one co- 
incide with the x-y plane and place the other at z = I. 
We approximate the light entering the cavity as a plane 
wave with phase velocity parallel to the z-axis. After 
each reflection inside the cavity, the reflected wave must 
have the same frequency as the incident wave, so we con- 
sider light of constant frequency p°. For simplicity, we 
set e — 1 in what follows. 

At leading order in the coefficients for Lorentz viola- 
tion, Eq. (|16| ) implies that the magnitude of the wave 
vector for each birefringent mode is \p± | = p± = [1 — 
(p ± <j)]p°. Decomposing the light entering the cavity 
into birefringent modes, we write 



E (x) = e- i P t [e i ^+ z (s n -E )i n 
+e^-*(£ T _.£ )e T 



(Dl) 



Here, p^± = [1 — (pj ± a^)]p a , where pj and a-] denote 
the values of p and a for light with phase velocity in 
the z direction. The unit vectors E\± are the associated 
birefringent basis. 

We suppose that the wave reflected at z = I has phase 
velocity in the —z direction. Decomposing this wave in 
the same fashion gives 



+e-^- z (i l ^E 1 )i^), 



(D2) 



where the subscript J, denotes quantities for phase veloc- 
ity in the — z direction. Similar expressions can be writ- 
ten for the electric field E n (x) after n reflections. For 
the nth reflection with n odd, the incident and reflected 
waves are related by E n+ i(x)\ z= i — e %Sl - n E n (x)\ z= i. A 
similar relation involving (5o, n holds for even n at z — 0. 



The complex factors e 



account for any phase 



change or loss due to transmission or absorption. They 
may depend on the interaction of the wave with the sur- 
faces and could also depend on the incident polarization 
and various coefficients for Lorentz violation. For sim- 
plicity, we suppose here that they are constant, and de- 
note them by Sq and Si. 

Superposing the contributions inside the cavity yields 
the total electric field as 

E(x) = e- vp0t [{e l ^+ z e ]+ e\ + + e l ^- z e^e\_) 

+e- lSo (e- ip t+ z i l+ i\ + + e- ip ^- z e^e\_)m] ■ M ■ E . 

(D3) 

At leading order, the matrix m is given by 



rn 



p + ercos£ crsin^ 
1 — i2p°l J ersin£ p — a cos £ 




(D4) 



where the bar signifies the average value over the two 
different propagation directions. The matrix M is the 
geometric series M = X)^Lo TO ™ = (1 ~ m)^ 1 . 

The resonant frequency is often viewed as the fre- 
quency at which a standing wave is produced in the cav- 
ity. However, this notion may fail for nonzero kp because 
the wavelength of light traveling in one direction can dif- 
fer from that for light traveling in the opposite direction. 
A more appropriate definition that applies also in the 
conventional case is to take the resonant frequency as 
the frequency maximizing the magnitude of the electric 
field or the energy density. The resonant frequency for 
a cavity is determined experimentally by measuring the 
transmitted light, so we adopt the energy density of the 
transmitted light as the relevant quantity. Using instead 



the magnitude of the transmitted electric or magnetic 
field yields the same result at leading order. 

We assume that the electric field Et of the transmit- 
ted light is proportional to the component of the total 
electric field in the cavity propagating in the z direction: 

E T cx e~ lp0t (e lp ^ z e n e\ + + e l ^- z e^e\_) ■ M ■ E . 

(D5) 

The time-averaged energy density is 

(u) = iRe (E* ■ D + B* ■ H) 

= \[E* ■{! + n DE ) ■ E + B* ■ (1 + n HB ) ■ B]. (D6) 

With this equation and the Faraday law ip°B = V x E, 
the energy density of the transmitted wave can be calcu- 
lated. Maximizing with respect to p and solving for p° 
yields the perturbed resonant frequency. We find 



+ a cos £ a sin £ ® \ E 
- = ^L- I JFiE? p-ctcos£ • -J- . (D7) 



5v E t 















En 



The barred quantities can be determined from 

asin£ = i(K _) n - \{k -) 22 - (k e+ ) 12 , 
ctcos£ = -\{k e+ ) lx + \{k e+ ) 22 - (k -) 12 , 



Kir 



(D8) 



which holds for a wave traveling in the +z direction. For 
a wave t rave ling in the — z direction, one must instead 
use Eq. (|d|) with sign changes for the parity-odd co- 
efficients: kdb —k-db, k-he —khe- The barred 



quantities in Eq. ( |32| ) then are merely those in (D8) with 
kdb = khe = 0. The net result is Eq. (p9), as desired. 



APPENDIX E: LABORATORY-FRAME 
QUANTITIES 



In the scenario of Sec. IV B 1 and in terms of the ma- 
trices k introduced in Eq. (|7|), the quantities ^0,1,2, 3, 4 
defined in Eq. @ can be written as 



A = i(l-3cos 2 x)[3e + (£ 



\ZZ 



-(«e-) 



ZZi 



-|(9e+ - e_)K t r 

+i/3 e [2(e_ - 3e+ cos 2 X ) sin q cos n @ T(k 0+ ) XY 
+ (3e+ — 2e_ — 3e+ cos 2 x) 
x (cos 77 cos n (B T(k 0+ ) xz + sin n e T(k 0+ ) YZ ) 
— 3e+(l — 3 cos 2 x) 

x (cos 77 cos n m T(k -) xz + sin fi e T(K _) yz )] 
+j€+(3 L sinx cos x(^o-) 



zz 



Ai = -i sinxcosx[3e + (K e+ ) 5 



-4, 



A* 



+fe+/3© sinx cos x[sin0 e T((K o _) KK - (k _) z ^) 

- cost] cos n e T ((k 0+ ) XY - {k _) XY ) 

+ smr 1 cosn (!) T((k 0+ ) xz - (£ -) XZ )] 

+±3 i [e_(* 0+ )^ 

+3e+(sin 2 x- cos 2 X )(k -) YZ ], 

-±sm X cos X [3e + {k e+ ) xz + e_(k e _) xz ] 

+fe+/3® sinx cos X [sinn (B T((k 0+ ) XY + (k _) XY ) 

+ cosi 1 cosn fB T((k -) xx - {k J) zz ) 

~sm7 1 cosn (B T((k 0+ ) YZ - (k -) YZ )} 

+\(3 L [e.{k 0+ ) xz 

+3e+(sin 2 x - cos 2 X ){k _) xz ], 

-\sin 2 X {3e + (k e+ ) XY + e_(k e _) XY ] 

-f e+/? e sin 2 X[sin n (S T{(k 0+ ) xz + (k _) xz ) 

+ cos?7 cos n B T({k 0+ ) YZ + (k q _) yz ) 

+ sin 7? cosfi e T((K _) xx - (k _) YY )} 

- |e+/3 L sin x cos X (k - ) XY , 



A 4 



■ X [3e + ({R e+ ) xx - (k e+ ) YY ) 



+e.((k e .) xx - (k e ^) YY )] 

+|e+/3 ffi sin 2 X [sin n B T((k 0+ ) YZ + (k _) YZ ) 

-cost] cos n m T({k 0+ ) xz + (k _) xz ) 

+2 sin 77 cos £L m T{k -) XY ] 



\e+l3 L sinx cos x[(«o- 



(«o-) 



(El) 



In this equation, we have introduced the quantities 



2 + e 
3e ' 



2-e 



(E2) 



which both reduce to 1 in the vacuum limit e — >• 1. 

The remaining coefficients are independent of e. The 
coefficients -Bo, 1,2, 3, 4 are given by 



Bo 
Bx 



Bo 



-\(3 L suv X (k 0+ ) XY ', 
-isinx[(^ + ) XZ -(« e -) A " Z ] 
+ i/3 e sinx[sinfi e T((^ 0+ ) xy + (ft„_) xy ) 
+ cos7 ? cos^ e T((K _) xx - («o-) ZZ ) 
-sin7 ? cosft ffi T((K 0+ ) rz - (K„-) yZ )] 
-i/3 L cosx((Ko+) XZ + (k„-) XZ ), 
isinx[(K e+ ) yZ -(A e _) yz ] 
-i/3 ffi sinx[sinft e T((£ _) yy - {k Q J) zz ) 
-cos?? cos n (B T({k 0+ ) XY - (k _) XY ) 
+ smr 1 cosn (B T({k 0+ ) xz - (k _) xz )} 
+y L cos X ((k 0+ ) YZ + (k _) YZ ), 



\XX 



B 3 = i C OSx[(Ke+) XX - (He-Y 

-(R e+ ) YY + (R^) YY } 
-±f3 S) cos X lsmn (B T((k 0+ ) YZ + (k _) YZ ) 



-cost/ cos n B T({k 0+ ) xz + (k„_) xz ) 
+2 sin 77 cos n (B T(k -) XY ] 
-±P L smx[(Ko-) XX -(«o-) yr ], 

S 4 = -|cOSx[(Ke+) Xr -(K e -)^] 

-±/?©cosx[sinft©T((£ 0+ ) xz + (k -) XZ ) 

+ cost] cos n B T{{k 0+ ) YZ + {h -) YZ ) 

+ sim 1 cosn m T((k _) xx - (k -) YY )} 

+ ±(3 L sm X (K -) XY . (E3) 

The coefficients Co, 1,2,3,4 are 

C = -|sin 2 X p e+ )^-(« e _)^] 

-\f3®sm 2 X [sintt (B T{{k 0+ ) YZ -Z{k _) YZ ) 
+ COST] cos {l(BT({k 0+ ) xz - i{k -) xz ) 
+2 sin 77 cos n (B T(k 0+ ) XY ] 
sinxcosx(K D -) zz , 

Ci = isin X cosxp e+ ) yz ~ («e-) yZ ] 

-i/3 e sin X cos X [sinr! ffi T((£ _) yy - {k _) zz ) 

-cost] cos n B T({k 0+ ) XY - {k _) XY ) 

+ sim 1 costt (S T{{R 0+ ) xz - (k _) xz )} 

+ lp L [(~ Ko+ ) YZ - (sin 2 X - cos 2 x){<-) Y % 

C 2 - isin X cosxp e+ ) xz - (^e-) XZ ] 

-i/? e sin X cosx[sinr!eT((Ko + ) Xy + {<-) XY ) 



C 3 



C 4 



- i<-) ZZ ) 

cos 2 x){<-) X % 
(r, \XYi 



+ cos 77 cos £l m T((k - ) xx 
— sin?]cosri0T((K o+ ) yz - 
+ i/3 L p 0+ ) xz -(sin 2 X - 

-|(l + cos 2 x)p e +) XF - ^e-7 j 
-i/? e (l + cos 2 x)[sintt ffi T((£ 0+ ) xz + (£ -) xz ) 
+ cos77Cosft©T((K 0+ ) yz + (k„_) 7Z ) 
+ sin 77 cos rJ e T((K _)- YX 
+ 1 /3 L sin x cos x(k _ ) XY , 
-i(l + cos 2 x)p e+ ) XX - 

-(K e+ )^ + (« e _)^] 

+ i/? ffi (l + cos 2 x)[sinf! e T((£ 0+ ) YZ + (« _) yz ) 
-cost? cos n (B T((k 0+ ) xz + (k _) xz ) 
+2 sin T) cos VL m T{k _ ) XY ] 

+ i/3 L sinxcosxpo-)* X - («o-) yy ]. (E4) 



APPENDIX F: SATELLITE-FRAME 
QUANTITIES 



- iK-) YY ) 

(Ke-) XX 



The quantities -A s ,o $s,c appearing in Eq. ( |50D of Sec. 
[V C 2 can be expressed in terms of the matrices k intro- 
duced in Eq. (Q) through the convenient combinations 



In terms of the various orientation angles specified 
in Appendix for the quantity A s we find 

A a = \ cos2C,[s\na{k e ,) xz - cosa{k e ,) YZ ] 
+\ sin2C[(l + sin 2 a)(k e ,) xx 
+ (1 + cos 2 a){k e ,) YY - sin 2a(k e ,) XY ] 
+ \P s [sina{{h ,) xz - (k ,) zx ) 

— cos a((k ,) Y z -(k Q ,) ZY )} 
+ j /3® [cos a cos 2( cos H®T 

x (cosr](k ,) YX — smr](k >) zx ) 

+ cosacos2Csinfi ffi T((K /) yy - {k ,) zz ) 

— sin a cos 2£ cos f2©T 

x(cosi 1 ((k ,) xx - (k ,) zz ) +smr,(k ,) ZY ) 

— smacos2£sinfi0T(K o /) J ' <:y 
+ cos 2 a sin 2£ cos n ffl T 

x (cos?7(k /) zx + smi](k i) YX ) 

+ cos 2 asin2Csinfi ffi T((K /) yz + (k ,) ZY ) 

+ sin 2 a sin 2£ cos f2©T 

x(cos7 ? ((w) XZ + (ko') ZX ) -sin77(K (/ ) Xy ) 
+ sin 2 asin2Csinfi e T(K ') zy 

— I sin 2a sin 2£ cos r2©T 

x (cos77(K o yz + sin77((K o XX - {<') YY )) 
-\ sin2asin2Csinfi ffi T(K /) xz ]. (Fl) 

The quantity A c is 

A c = —4 cos£sina(/t e ') yz 

— j cosCcosa(«; e ') xz 

- 1 sinCsin2a((K e /) xx - (k e ,) YY ) 
+i sinCcos2a(K e /) xy 

+ 1 /3 s [sinC((^) Xy -(^) YX ) 
+ cosCcosa((K /) zx - {k ,) xz ) 

+ cosCsina((^) Zy - (^) 7Z )] 

+ 1 /3© [cos a cos C cos r2©T 

x(cos?7((k o xx - (^o') ZZ ) +shi77(K o zy ) 

+ cos a cos C sin n®T(ft < ) xy 

+ sin a cos C cos Q®T 

x (cos 77(k - ) YX — sin 77(^0' 

+ sinacosCsinfi e T((K /) yy - (k ,) zz ) 

+ cos 2a sin £ cos f2©T 

x(cos77(K o yz + sin77((K o XX - (k ') YY )) 
+ cos 2a sin ( sin fl e T(k ol ) xz 

— sin 2a sin £ cos f2® T 



sini 



x (cos?7(k /) 
+ sin 2a sin £ sin f20T(K o < ) 



r;((^) XY + (^) yx )) 



YZ1 



(F2) 



The quantity B s is 



B s 



■ sin£sina(K e ') 



YZ 



| sin C cos a(k e i ) xz 



xx 



16 



{He') 



YY\ 



i£sin 2a((«v) 
+ | cos C cos 2a(k e ,) XY 
-|/3 s [sinCcosa((K o XZ + (k <) zx ) 
+ sinCsina((K o rZ + (k >) ZY ) 
-cosCsin2a((K ') XX - {^o') YY ) 
+ cosCcos2a((K o Xy + (fto>) YX )] 
+|/?e [— cos a sin^cosfi®T 
x(cos7?((/t ,) xx - (k >) ZZ ) +smr](k ,) ZY ) 

— cos a sin £ sinQ^T (iio>) XY 

— sin a sin £ cos fl^T 

x (cos?7(K ') rx — sim](k i) zx ) 
-sinasinCsinft ffi T((K /) yy - (k ,) zz ) 
+ cos 2a cos C cos f2©T 

x (cosr7(K o rz + sin7?((/v)* x - (« ') yy )) 
+ cos 2a cos C sin f2©T(fiv ) xz 

— sin 2a cos £ cos f2©T 

x (cost7(k o XZ - sm V ((k ,) XY + {k ,) YX )) 



+ sin 2a cos C, sin f2©T(K ' ) YZ ] 
The quantity 2? c is 



B c = ^(cos 2 a- sin 2 a cos 2 ( + sin 2 C)(k 6 /) xx 



(F3) 



+ (sin 2 a — cos 2 a cos 2 £ + sin 2 ()(k e i) YY 
-^sin2a(l + cos 2 ()(k e ,y 



\XY 



16 1 

3 



+ ^sin2C[sina(K e /) AZ - cosa(fv) ] 
-^/3 s [2(cos 2 a - sin 2 a cos 2 C + sin 2 Q{k ol ) xx 
+2(sin 2 a - cos 2 a cos 2 ( + sin 2 C)(k > ) YY 
+ sin 2a(2 cos 2 C + sin 2 ()((k ,) XY + (k ,) YX ) 
+ smasm2(((k ,) xz + (k a ,) zx ) 
-cosasm2(((k ,) YZ + 



+ 



Ki * 



, [sin 2a(l + cos 2 Q cos fi©T 



x^osiji^f 2 + smj 1 ({k ,) xx - (k a ,) YY )) 
+ sin2a(l + cos 2 Q sin fl m T(k ,) xz 

— sin a sin 2£ cos Q&T 

x(cos V ({k ,) xx - {k ,) zz ) +sin V (k ,) zx ) 

— sinasin2£sinf2©T(K /) xy 
+ cos a sin 2Q cos fi©T 

x (cos?7(K /) yx - smrj(k r) zx ) 

+ cosasin2£sinf2©T((K /) yy — (k >) zz ) 

— 2(sin 2 a cos 2 C — cos 2 a) cos fi©T 

x (cosr](k r) xz - smrj(k r) XY ) 



— 2(cos 2 a cos 2 ( — sin 2 a) cosf2©Tsin77(K ') yx 
— 2(cos 2 a cos 2 C — sin 2 a) sinf2©T(K /) yz 
+2 sin 2 C(cos 0©T cos rj{k ,) zx 

+ s\nn <B T{k ol ) ZY )}. (F4) 

Finally, the quantity C is 



C = yg(3sin 2 a sin 2 £ — l)(K e ') 
+ ^(3cos 2 asin 2 C- l)(k e >) 



yy 



xy 



+ ^(3cos 2 C- 1)(k c /) zz + ^|sin2asin 2 C('v) 

sin2C[sina(K e /) xz - cosa(k e >) YZ ] 
+ |/3 S [(3 sin 2 a sin 2 ( - l)(k a ,) xx 
+ (3 cos 2 a sin 2 (- l)(k >) YY 
+ (3cos 2 C-l)(«oO ZZ 
-|sin2asin 2 C((^) Xy + («o') yX ) 
+| sin2C(sina((K o XZ + {k >) ZX ) 
-cosa((ko') YZ + (Ko') ZY ))] 
+ |/3©[(3sin 2 asin 2 C- l)cosO ffi T 
x (cosri(k ,) xz - sin7y(K /) xy ) 
+ (3 cos 2 a sin 2 ( - 1) (cos 0©T sin r](k ,) YX 
+ sinn e T(k ,) YZ ) 
-(3cos 2 C - l)(cosfi©Tcosry(K -) zx 
+ sinf7 ffi T(K o zy ) 
— | sin 2a sin 2 £cosf2©T(cosr7(«; ') yz 

+ sinr;((^) XX -(^) yy )) 
-| sin 2a sin 2 C sin n m T(k <) xz 
-|sinasin2Ccosfi ffi T(cos?7((K /) xx - (k >) zz ) 
+ smr](k ,) ZY ) 

-§ sin a sin 2( sin fi©T(K /) xy 

+| cosasin2(cosri0T(cos77(K o ') yx 

-sinr7(K o ZX ) 

+ |cosasin2Csinfi ffi T((K o yy - (k & ) zz )\. (F5) 
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